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Accurate, efficient, automatic methods for computing the complex error function to any precision are 
detailed and implemented in an American Standard FORTRAN subroutine. A six significant figure table of erfc 
2, e' erfc z, and e' erfc( — z) is included for z in polar coordinate form with the modulus of z ranging from to 
9. The argand diagram is given for erf z. 
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1. Introduction 

In computing many of the functions of mathematical physics, for example, Fresnel integrals, Dawson's 
integral, Voigt function, plasma dispersion function, etc., difficulties are frequently encountered. Since 
these functions may be expressed in terms of the error function of complex argument, we have chosen this 
function for Part IV. 1 The major part of the coding of the power series, continued fraction and asymptotic 
expansion computations for complex arguments will carry over equally well for other functions. 

As Part I was devoted to the error function of a real variable, the probability function and other related 
functions, Part IV will only emphasize those functions and pitfalls due to complex arguments. 

While accuracy over the entire domain of definition remains our main concern, the methods employed 
ensure efficiency, portability and ease of programming and modification. 

If one supplies approximate values for the maximum machine value, minimum machine value, the upper 
bound of the sine, cosine routine, and the upper bound to the acceptable relative error and gives the square 
root of 7r to the required number of significant figures, the detailed methods will work for computations 
ranging from very low precision to multi-precision. 

The argand diagram of erf z is included as well as the implementing ANS FORTRAN program and a six 
significant figure table of erfc z, e* erfc z and e* erfc( — z) for z in polar coordinate form with the modulus 
of z ranging from to 9. 

2. Mathematical properties 

Relevant formulas are collected here for completeness and ease of reference. In keeping with the conven- 
tion of the Handbook [l], 2 z = x + i y is a complex variable. 



•Mathematical Analysis Division, Center for Applied Mathematics. 

1 Part I. Error, Probability, and Related Functions. J. Res. Nat. Bur. Stand. (U.S.). 74(3): 211-224; 1970, Part II. The Exponential Integral £.(*). J. Res. 
Nat. Bur. Stand. (U.S.). 78(4): 199-216; 1974. Part III. The Sine, Cosine, Exponential Integrals, and Related Functions. J. Res. Nat. Bur. Stand. (U.S.). 
80(2): 291-311; 1976. 



* Figures in brackets indicate literature references at the end of this paper. 
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A. Definitions 

erf z = — -, — I * e~' dt 
Vtt ° 

erfc z = — -, — | °° e~' dt = 1 — erf z 

yJTT ' 

(The path of integration is subject to the restriction arg / — a with |a| < -r- as t — oo along the path. If 
Rt 2 remains bounded to the left, a = -r- is permissible.) 

w(z) = e'* 1 (1 + — -, — \ z e' 2 dt) = e~* 2 erfc ( -iz) = e* 2 erfc £ (£ = ~**) 
i fa. e~* efa 2i2 f o. e~* dt , r .. 

= ir!-^r = — '»i^f ( ' i>0) 

with F(z) = e~* 2 J* e' 2 dt (Dawson's Function) 

i -< 2 j 

and Z(z) = — -. — \°° (Plasma Dispersion Function) 

yir t— z 

E(z) = C(z) + i 5(«) = j; e ™" dt = \l cos (-^-) dt + i \\ sin {-^-) dt 
= ~2~ erf [— g - (1 ~0 z] (Fresnel Functions) 

F fe'> = ?T7wi II £ TT^ = (-^7-) ,/2 ^ «fc u 

(47T/)" 2 °° U 2 +0tf 2 4/ 



where 



oj= (l-n:)/2 * 1 ' 2 
= U(x,t) + i V[x,t) (Voigt Function) 

1 ( . e- (Jf - y)2/4t i yg-c-y) 2 ^! ^ 



(4rt) iy!j - 1+y ^ (4ttO i/2 J -°° l+y 
e"' dt 



it/ . a foo c~' af 1 

lUcu; - — ].„, ( u _ f j» + fl2 = — J^ U(u/a, l/4a 2 ) (Line Broadening Function) 



B. Series Expansions 

erf z = -?- £ ( ~ / 1)V '' 41 
Vtt «=° n! (2n + 1) 

_ 2 2 - 2"z 2 " +1 

e I, 



V?r -o 1-3. ..(2n + l) 

C. Continued Fraction (Rz>0) 
e' 2 erfez= * U \!L _i_ _KL _2_ ] 

V7T L 2 + 2+ 2+ 2+ z + ' "J 

= 2z \ 1 1»2 3*4 -. 

Vtt l 2**+l - 2^ + 5 - 2z 2 +9 -•••■! (" Even "Form) 
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D. Asymptotic Expansion 

e erfc2 ~^r L T + £ i^f 1( ' ' g ' < "T" ) 

E. Symmetry Relations 

erf (-2) = - erf 2 



erf 2= erf 



w(~ z)= 2e~* — w(z) 



w(z)= w{-z) 
C(-z)= - C(z),S(-z)= -S(z) 
C(iz)= i C(z), S(iz) = —i S(z) 

C(7) = C(7),S(7) = m 

3. Method 

The main functions under consideration are the error function ERFZ, the complementary error function 
ERFCZ, and the exponential of z 2 times the complementary error function EZ2CZ. All other functions may 
be obtained from these three. To simplify testing, computations are performed for z in the first quadrant AZ 
and symmetry relations are then employed to make adjustments for other quadrants. For the special case z 
= 0, no computations are performed and the following function values are returned: ERFZ = 0, ERFCZ = 
1 and EZ2CZ = 1. 

Real type variables are used throughout to readily allow for double precision computation if greater accu- 
racy is needed. The machine dependent constants are placed in a labeled section at the beginning of the 
subroutine. Function references are likewise grouped together when possible and attention called to the 
statement labels of the remaining function references. Real and imaginary parts of complex variables have 
R and I as final characters. 

Since EZ2CZ for z in the first quadrant is machine representable even with the real and imaginary parts 
of z equal to the maximum machine value CMAX (provided its reciprocal is larger than the minimum 
machine value CM1N), checking for the range of the argument z has been omitted. However, the extensive 
range necessitates a fair amount of testing for overflows. Underflows are assumed to be set to zero. 
Overflows are set equal to the maximum machine value and an error indicator IERR set for the number of 
functions affected. If only EZ2CZ lies outside the machine range, IERR = 1, otherwise IERR = 3. As often 
as possible, computations are arranged so as to give the correct results for the three functions if they lie 
within the range of the machine. 

In computing the modulus RHO = V(AZR) 2 + (AZI) 2 of a complex quantity AZ = AZR + i AZI in the 
first quadrant, RHO may lie in the machine range but (AZR) 2 , (AZI) 2 or their sum may be outside the range. 
We select the larger ARIMX and smaller ARIMN of either AZR or AZI, and compute the ratio RMNMX = 
ARIMN/ARIMX. A factor of RHO called PRHO is computed as the square root of (RMNMX'RMNMX + 
ONE). This factor, which is greater than or equal to one and less than or equal to the square root of 2, can 
then be used to check for overflow. The quantity ARIMX must be less than CMAX divided by PRHO for 
RHO to lie in the machine range. A similar procedure is followed in computing the real and imaginary parts 
of (AZ) 2 = Z2R + i AZ2I = (AZR) 2 - (AZI) 2 + i 2AZR*AZI with first checking to ensure ARIMX is 
greater than or equal to 1. 

Analysis has indicated and testing confirmed that the power series PS is most useful from the standpoint 
of accuracy and efficiency for RHO less than RHOLS( = 1.5) and when AZR is less than or equal to 1 provid- 
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ed RHO is less than AELL (=V - ln(TOLER)) where TOLER is the upper limit for the relative error. The 
continued fraction expansion CF is most useful for AZR greater than 1 and RHO greater than or equal to 
RHOLS. The asymptotic expansion AE is most useful for AZR less than or equal to 1 for RHO greater than 
or equal to AELL. For RHO greater than or equal to RHOLC (=V0.5/TOLER) a rearrangement of arith- 
metic operations for the first term of the asymptotic expansion is necessary to maintain the accuracy of 
EZ2CZ. In the continued fraction and asymptotic expansion regions only, EZ2CZ is first computed; it tends 
to zero for large \z\ and the exponential of — z 2 tends to infinity for small AZR, To maintain accuracy 
here, we compute the exponential of — Z2R/6 and do continuous multiplication and testing with appropriate 
factors to obtain ERFC. The imaginary part of (AZ) 2 is tested against ULSC, the upper limit of the sine, 
cosine routine. 

Figure 1 below maps the regions for the various methods. 

v. - ■ *-C 



6 _ 



5 - 

AELL 
h _ 



3 - 
2 - 
1 - 



AE 



PS 



1 



CF 



Figure 1. Parameter plane. 

The dividing line RHO = AELL between the use of the power series and the asymptotic expansion and 
RHO = RHOLC are the only boundaries subject to the required precision. Single and double precision 
results, for example on the Univac 1108, are the results of two different methods in the region where RHO 
roughly lies between 4.3 and 6.4. 

This mapping of the region ensures for the required precision that the least number of terms are com- 
puted and the loss of significance is kept to a minimum. While the second form of the power series is prefer- 
able for real positive z, since all terms are positive, the first form ensures greater accuracy for complex z 
since the real and imaginary parts of the terms may be positive, negative or zero for any RN( = n). The power 
series is there more rapidly convergent; the relative error may then be approximated by the ratio of the term 
to the sum of terms; comparison of this approximant with the tolerance for the relative error controls the 
number of terms needed. Since the terms tend to zero through underflow, there will always be convergence 
even if the tolerance is made considerably smaller than the precision of the machine. The power series is 
evaluated using recurrence relations in the following form: 

TWO RNF ™n 

ERFZ = XrH^L.^ £ SGN(RN)*TM(RN) = ^^~ * SUM 



SQRT(PI) RN - =0 
where 

SGN(O) = 1, SGN(RN + 1) = -SGN(RN) 
TM(RN) = ((AZ**(2*RN + 1))/1*2* . . . RN)/(2*RN + 1) 
TM(RN) = PTM(RN)/DN(RN) 

PTM(O) = AZ, PTM(RN + 1) = (AZ**2)*PTM(RN)/(RN + 1) 
DN(0) = 1, DN(RN + 1) = DN(RN) + 2 



RTPI 
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In determining the terminal value of RN, normalization is necessary to avoid overflows and underflows. 
The normalization factor TMAX is the maximum of the absolute value of the real and imaginary parts of 
TM and SUM. If TMAX equals zero or |TM/TMAX| 2 underflows, RNF = RN. 

If |SUM/TMAX| 2 underflows, additional terms are obtained. Otherwise, if | TM/TMAX 1 2 I | SUM/ 
TMAX | 2 is less than TOLER 2 , then RNF = RN. 

The continued fraction expansion starts to converge more slowly as z tends to zero. The "even" form is 
used since the required number of terms is halved at the expense of very little extra computation for suc- 
cessive numerators and denominators. The continued fraction is evaluated by using the recurrence relations 
in the "forward" direction. The number of terms needed is determined by checking to see if the relative 
error of two successive convergents is less than the tolerance. On the other hand, if the relative error 
remains constant or starts to increase, the recurrence is terminated and the prior convergent taken as the 
value of the continued fraction. In this way, the process is always terminated when maximum precision is 
attained. 

The "even" form of the continued fraction takes on the following implementation: 

RNF 
EZ2C Z=4^-RN n =l AM(RN) 



RTPI BM(RN) 

with AM(1) = 1, AM(RN + 1) = -WM(RN + 1)*(WM(RN + 1) + 1) 

BM(1) = 2*(AZ**2) + 1, BM(RN + 1) = BM(RN) + 4 
where WM(1) = - 1, WM(RN + 1) = WM(RN) + 2 

EZ2CZ = (AZ*(FM/GM))*2/RTPI = (AZ*F(RN))*2/RTPI 
where FM(-l) = 1, FM(0) = 

GM(-l) = 0, GM(0) = 1 
and FM(RN) = BM(RN)*FM(RN - 1) + AM(RN)*FM(RN - 2) 

GM(RN) = BM(RN)*GM(RN - 1) + AM(RN)*GM(RN - 2) 

The relative error may be approximated by [F(RN) - F(RN - 1)]/F(RN) = RE(RN). If the modulus 
squared of the relative error REM2(RN) is less than the square of the tolerance divided by 8, RNF = RN. If 
REM2(RN) is greater than or equal to REM2(RN - 1), then RNF = RN - 1. Normalization is likewise 
necessary here to avoid overflows in computing the relative error and its modulus squared and also in the 
generation of the successive convergents. 

The asymptotic expansion is likewise evaluated using recurrence relations in the following form: 

i RNF 

EZ2CZ = c *, m , £ SGN(RN)*TM(RN) = SUM/RTPI 
byRT(PI) rn=0 

where SGN(O) = 1, SGN(RN + 1) = -SGN(RN) 

TM(RN) = (1/AZ)*(P3* . . . (2*RN - 1))/(2*(AZ**2))**RN 

TM(0) = 1/AZ, TM(RN + 1) = DN(RN + 1)* TM(RN)*(1/(2*(AZ**2))) 

with DN(1) = 1, DN(RN + 1) = DN(RN) + 2. 
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The relative error may be approximated here by the ratio TM/SUM. The convergence test precedes the 
divergence test and is implemented as REM2 less than (TOLER**2)/8 to attain greater accuracy in both the 
real and imaginary parts. If the modulus squared of the term remains the same or increases, the prior sum is 
taken as the final sum. 

For z along the imaginary axis, the error function is purely imaginary; the real part of erfc z = 1 and of 
e* 2 erfc z — e" (AZ " 2 . No difficulties arise in the use of the power series. However, since the asymptotic expan- 
sion is given for \z\ — oo, the correction must be applied for AZR — 0. 

The following table gives an indication of the number of terms needed to obtain maximum machine accu- 
racy on the Univac 1108 with the various methods of computation. 



Method 



Number of Terms 



Single Precision 
TOLER = .745E-8 



Double Precision 
TOLER = .867D-18 



Power Series 
Continued Fraction 
Asymptotic Expansion 



50 
25 
22 



112 
99 

45 



4. Range 

If the real part of z is zero or positive, e" 2 erfc z is valid for z throughout the entire machine range. Other- 
wise, the real part of z 2 is essentially limited by the range of the exponential library subroutine with the 
imaginary part of z 2 limited by the range of the sine, cosine library subroutine. 

5. Accuracy and Precision 

The maximum relative error, generally in erfc z, except for regions in the immediate neighborhood of 
zeros of the real and/or imaginary parts of the functions is 8E-6 for single precision computation on the 
Univac 1108. 

The precision may be varied by changing the value of TOLER. 

6. Timing (Univac 1108 Time/Sharing Executive System) 

The time estimates given below are highly dependent on the operating system environment and conse- 
quently should not be relied on for critical timing measurements. 





Time (Seconds) 






Single Precision 


Double Precision 


Region 


TOLER = .745E-8 


TOLER = .867D-18 


ZR = 0(.1H, ZI = 0(.2)8 (1681 values) 


5.94 


22.5 


Method 


Maximum Time/Evaluation 


Power Series 


.0101 


.038 


Continued Fraction 


.0088 


.052 


Asymptotic Expansion 


.0035 


.0093 



7, Testing 

The language of the subroutine was checked for conformity with the PFORT VERIFIER. 3 Test argu- 
ments were devised and used in the analysis of the subroutine with the PROFILER. 4 



'The PFORT Verifier, A, D, Hall and B. G. Ryder, Bell Laboratories, Murray Hill, NJ. Proceedings of the Computer Science and Statistics Eighth 
Annual Symposium on the Inlerfaee, University of California, Los Angeles, February 13-14, 1975, 

♦ Program Execution Profiles, C. Sande, World Bank, Washington, D.C. Proceedings of the Computer Science and Statistics Eighth Annual Symposium 
on the Interface, University of California, Los Angeles, February 13-14, 1975. 



666 



The subroutine was used to obtain related functions which were checked against available published 
tables ([l]-[5], [7]-[9], [13], [15]). Single precision results covering the 9 X 9 grid were compared against 
double precision results. This precision test particularly verified that the scaling operations were valid and 
undetected overflows had not occurred. 

Where applicable, results were obtained by independent methods, for example, the power series and con- 
tinued fraction, the power series and asymptotic expansion or all three methods. Otherwise, both forms of 
the power series or continued fraction were used. In addition, numerical integration of various forms of the 
integral representation was employed. 



8. Test Values 

Six-significant figure tables of erfc z, e* 2 erfc z and e' 2 erfc( — z) are included in the appendix. The com- 
putations were performed with double precision arithmetic to provide more accurate values for checking 
purposes. The arguments are in polar coordinate form z = Qe" for q = 0(.02).2(.1)3(.5)9, 9 = 0°(15°)30°, 
37.5°, 45°, 50°(10 o )90°. Values of the functions for z in other quadrants are readily obtainable with sym- 
metry relations. 



9. Argand Diagram of erf z 



_ 2 

V7T 



ir 



dt 



Since erf z and erfc z are complementary functions, erf z is not tabulated but the argand diagram of the 
function is included to illustrate divergent and convergent spirals. 



erf ~, 8=0° 

Convergent spirals, 8<45° 

Cornu's spiral of Fresnel Integrals, 8=45,' 

Divergent spirals, 6>A5° 

3£ / p et 2 dt, 6=90° 



\ 



-It 




\, 



. H.T. 



*&i 



^ ^-ir" / 1 I 



f 



/ 



P. erf z 
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10. Special Values 



Relevant values are collected here for completeness and ease of modification and checking of the 
program. 



Zeros 



erf z f =0 



z, 



1.45061616 + i 1.88094300 
2.24465927 + i 2.61657514 
2.83974105 + i 3.17562810 



s z, 

1 1.99146684 - i 1.35481013 

2 2.69114902 - i 2.17704491 

3 3.23533087 - i 2.78438761 



C(z,)=0 

*. 

1.7437 + i .3057 
2.6515 + i .2529 
3.3204 + i .2240 



S(z,)=0 

z, 

2.0093 + i .2885 

2.8335 + i .2443 

3.4675 + i .2185 



Maxima and Minima of Fresnel Integrals 



s M, = aV47Tl) 

.779893 

1 .640807 

2 .605721 

3 .588128 



m, - C(V47+3) 

.321056 

.380391 

.404261 

.417922 



M* s = S(V4TF2) 

.713972 

.628940 

.600362 

.584942 



m* = 5( x /47+4) 

.343416 

.387969 

.408301 

.420516 



Radius of Univalence q 
q = 1.5748376 .92413887 

Maximum and Inflection Point for Dawson's Integral 
F(.92413 88730) = .54104 42246 
F(l. 50197 52682) = .42768 66160 
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Related Constants 



1° 


= 


1.74532 


92519 


94329 


57692 


36907 


68488 


61271 (-2)r 


7T 


= 


3.14159 


26535 


89793 


23846 


26433 


83279 


50288 


tt/2 


= 


1.57079 


63267 


94896 


61923 


13216 


91639 


75144 


Vtt 


= 


1.77245 


38509 


05516 


02729 


81674 


83341 


14518 


V*/2 


= 


.88622 


69254 


52758 


01364 


90837 


41670 


57259 


2tt 


= 


6.28318 


53071 


79586 


47692 


52867 


66559 


00576 


2Vtt 


= 


3.54490 


77018 


11032 


05459 


63349 


66682 


29036 


e 


= 


2.71828 


18284 


59045 


23536 


02874 


71352 


66249 


1/tt 


= 


.31830 


98861 


83790 


67153 


77675 


26745 


02872 


2/tt 


= 


.63661 


97723 


67581 


34307 


55350 


53490 


05744 


1/Vtt 


= 


.56418 


95835 


47756 


28694 


80794 


51560 


77258 


2/Vtt 


= 


1.12837 


91670 


95512 


57389 


61589 


03121 


54517 


l/2ir 


= 


.15915 


49430 


91895 


33576 


88837 


63372 


51436 


1/2 Vtt 


= 


.28209 


47917 


73878 


14347 


40397 


25780 


38629 


1/e 


= 


.36787 


94411 


71442 


32159 


55237 


70161 


46086 


erf 1 


= 


.84270 


07929 


49714 


86934 


12206 


35082 


60926 



\y* = (Vx/2) erf 1 = | -^^ 



(VTr/^erf 1 = E 



2" 



o 1.3. ..(2n + l) 



' * 2 



1 



lo e ' dt - So nl(2n + \) 

t l 

«>o (2n)!(4n + l) 

1 
& (2n + l)!(4n + 3) 

■» 1-3. . .(4n + l) 

So 1.3. ..(4n + 3) 

e U e at - ko 1.3. ..(2n + l) 



= .74682 41328 12427 02539 94674 36131 85300 



= 2.03007 84692 78704 97553 90899 25665 95044 



= 1.46265 17459 07181 60880 40485 86856 98815 



= 1.10473 79393 59804 31710 17580 11494 42058 



= .35791 38065 47377 29170 22905 75362 56757 



= 1.28407 89880 95736 69733 77386 73036 75360 



= .74599 94811 82968 27820 13512 52629 19684 



= .53807 95069 12768 41913 63874 20407 55675 



Typical Tolerances and Their Natural Logarithms 



2 -I6 
2 -24 
2 -J7 
2 -36 
2 -<»8 
2 -56 
2 -60 
2-108 



0.15258 78906 25(-4) 

.59604 64477 53906 25(-7) 

.74505 80596 92382 8125(-8) 

.14551 91522 83668 51806 64062 5(-10) 

.35527 13678 80050 09293 55621 33789 0625(-14) 

.13877 78780 78144 56755 29539 58511 35253 90625(-16) 

.86736 17379 88403 54720 59622 40695 95336 91406 25(-18) 

30814 87911 01957 73648 89564 70813 58837 09660 96263 71446 21112 

38390 20729 06494 14062 5(-32) 



log.(2- 16 ) = -11.09035 48889 59124 95067 57139 43330 82508 
log,(2-») = -16.63553 23334 38687 42601 35709 14996 23763 
log.(2-") = -18.71497 38751 18523 35426 52672 79370 76733 



669 



log.(2- 36 ) 
log,(2- 48 ) 
log.(2- 5 *) 
log.(2- 60 ) 
log«(2- ,os ) 



-24.95329 85001 58031 13902 03563 72494 35645 

-33.27106 46668 77374 85202 71418 29992 47526 

-38.81624 21113 56937 32736 49988 01657 88781 

-41.58883 08335 96718 56503 39272 87490 59408 

-74.85989 55004 74093 41706 10691 17483 06935 



Maximum and Minimum Machine Values and Their Natural Logarithms 
(Univac 1108 Single and Double Precision Limits) 



NBC -Number of binary digits in the (biased) characteristic of a floating point numb 



er 



2 127 

2-129 

log,(2»") 
log*(2- J29 ) 



NBC = 8 

0.17014 11834 60469 23173 16873 03715 88410(39) 
.14693 67938 52785 93849 60920 67152 78070(-38) 
88.02969 19311 13054 29598 84794 25188 42414 
-89.41598 62922 32944 91482 29436 68104 77728 



2»023 
2">025 

log.(2 10 ") 
log,(2- ,02s ) 



NBC = 11 

0.89884 65674 31157 95386 46525 95394 51236(308) 

.27813 42323 13400 17288 62790 89666 55050(-308) 

709.08956 57128 24051 53382 84602 51714 62914 

-710.47586 00739 43942 15266 29244 94630 98227 
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C**** APPENDIX 

IMPLEMENTING PROGRAM 
C LANGUAGE. AMERICAN NATICNAL STANDARD FORTRAN 
(SUBROUTINE SUBJECTED TO PFORT VERIFIER*) 
DEFINITIONS. Z, A COMPLEX V AR I ABLE = ZR-M ZI 

ERF(Z) = (2/S0RT(PI) )* INTEGRAL( EXP( -T**2 > )DT FRCM TO Z 

=ERFZfi+I ERFZI 
ERFC<Z>=(2/SQPT(FI> > *INTEGRAL ( EXP { -T**2 ) > DT FROM Z TO 

INFINITY 
=1-ERF(Z> 
=ERFCZfi+l ERFCZI 
EXP(Z**2)*ERFC< 2>=EZ2CZR+I EZ2CZI 
SYMMETRY RELATIONS 
ERF<-Z)=~ERF<Z> 
ERFCZ CCNJG>=CCNJG(ERF(Z) ) 
ERFC<-ZJ=2-EfiFC(Z) 
ERFC(Z CCNJG1=CCNJG(ERFC(Z) » 

EXP(Z**2)*ERFC(-Z>=2*EXP(Z**2)-EXP(Z**2>*ERFC(Z) 
EXP(Z CCNJG**2>*ERFC(Z CON JG) =CCN JG( EXP( Z**2 » *ERFC< Z > ) 
SPECIAL CASEi Z=0 
ERF<2>=0 
ERFC<Z>=1 

EXP<Z**2)#ERFC(Z)=1 
USAGE. CALL ERRZ < ZR, ZI , ERFZR .ERFZ I, ERFCZR, ERFCZI .EZ2CZR. 

EZ2CZIt IERR > 
ARGUMENTS 

(REAL TYPE VARIABLES ARE USED THROUGHOUT TC READILY 
ALLCh FOR DCUELE PRECISION COMPUTATION. REAL AND 
IMAGINARY PARTS CF COMPLEX VARIABLES HAVE R AND 
I AS FINAL CI-ARACTERS.) 

ZR.ZI REAL(OR DOUBLE PRECISION) TYPE INPUT 

ERFZRtERFZI (SAME TY"E AS Z) OUTPUT 

ERFCZR.EFFCZI •» OUTPUT 

EZ2CZR.EZ2CZI " OUTPUT 

IERR INTEGER TYPE OUTPUT 

IERR 

NORMAL RETURN 

1 EXP(Z**2)*ERFC(Z) INVALID 

2 ERF(Z) .ERFC(Z) INVALID 

3 ERF(Z)iERFC(Z) tEXP ( Z** 2 J *ERFC < Z > INVALID 
(Z IN 2ND OR 3RD QUADRANTS(ZR .LT. 0)) 

COMMONLY ISED INTERNAL VARIABLES 

AELL LOWER LIMIT OF tZl FOP ASYMPTOTIC 

EXPANSION A.E.) ABS(ZR) .LE. 1 
AZI AbSCZIJ 
AZR A6S(ZR> 
AZ2I AES( IMAG(Z**2> 1 
CMAX MAXIMUM MACHINE V*LUE 
CMIN MINIMUM MACHINE VALUE 
PEM2 MODULUS SQUARED OF RELATIVE 

ERROR(R.E. ) 
1R.E. 1**2=REM2(N ) 
REPM2 REM2(N-1> 
RHO !Z1=S0RT(ZR**2+ZI**2) 
RHCLC LOWER LIMIT OF RHO FOR USING ONLY 

1ST TERM Op A.E. 
KHOLS UFPER LIMIT OF RHO, UNRESTRICTED 
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56 C AeS(ZR). FGR POWER SERIES 

ES C RTPI S0RT{PI=3.14... ) 

60 C SUMM2 !SUMI**2=SUMR**2+SUMI**2 
£1 C OR !SUM/TMAXI**2 

tZ C TMAX NORMALIZATION FACTOR 

£3 C TMM2 !TM!**2=TMR**2+TMI**2 

£4 C OR !TM/TMAX!**2 

C5 C TOLER UPPER LIMIT FOR RELATIVE ERRORS 

€6 C TCLEF2 TGLER**2 

t7 C T0L2 TCLER2/8 

ee C ULSC MAXIMUM ARGUMENT FOR SIN/COS ROUTINE 

f9 C MOCIFICATICNS. 

70 C THE CODE IS SET UP FOR SINGLE PRECISION COMPUTATION 

71 C ViITH SINGLE PRECISION FUNCTION REFERENCES AND SINGLE 

72 C PRECISION MACHINE DEPENDENT CONSTANTS. FOR THE UNIVAC 

73 C 1108t CNAX APPFCX. 2**1 27, CMIN= 2** ( -* 1 29 ) »ULSC=2**20 AND 

74 C TOLER=.745E-e. FTPI IS GIVEN IN DOUBLE 

75 C PRECISION FORMAT TO 19 SIGNIFICANT FIGURES. 

76 C DOUBLE PRECISION RESULTS ARE OBTAINED BY INSERTING 

77 C (1) THE DCUBLE PRECISION TYPE STATEMENT 

78 C (2) DCUELE PRECISION INTRINSIC FUNCTION REFERENCES - 

79 C DABStCMAXl ANC CMIN1 

80 C (3) DOUBLE PRECISION EXTERNAL FUNCTION REFERENCES - 

61 C DCCStDEXP.CLOGtDSIN AND DSQRT AND 

£2 C (4) FOR ThE UNIVAC 1108 ADJUSTING THE CONSTANTS 

63 C CMAX APPPOX. 2**1023 , CM IN=2** <-l 025 ) ,ULSC=2**56 AND 

64 C TCLEfi=.ee7C-18. 

65 C THE DETAILEO NETFCDS SHOULD WORK FOR ANY PRECISION 

66 C IF THE MAChlNE CEFENDENT CONSTANTS ARE CHANGED 

67 C WITH RTPI GIVEN TC THE REQUIRED NUMBER OF SIGNIFICANT 
ee C FIGURES. 

89 C METHOD. 2=ZR+I ZI = RHC*EXP < I * AR CTAN ( Z I/ZR > ) 

90 C ALL METHODS APPLY TO AZ=ABS(ZR)+J AES(ZI)=AZR+ 

91 C I AZI. USE IS THEN MADE OF SYMMETRY RELATIONS. 

92 C POWER SERIES 

93 C RHC .LT. RH0LS(=1.5) 

94 C AZR .LE. 1, RHOLS .LE. RHO .LT. AELL 

95 C AELL=SQRT(-LOG(TOLER) ) 

96 C ERF(AZ)=(2/S0RT(PI) )*SUM(SGN(RN)*TM(RNI ) 

97 C PN=0.1«...,RNF 

98 C SGN(0I=1 

59 C «GN(RN+1 )=-SGN(RN) 

ICO C TM(RN) = ( ( AZ**(2*RN+1 ) )/l*2...RN)/(2*RN+l ) 

101 C TM(RN)=PTM(RN)/DN<RN) 

102 C PTM(0)=AZ 

103 C FTM(RN+1 )=( AZ**2 ) *PTM < RN > / < RN+1 ) 

104 C DN<0)=1 

105 C DN(RN+1 )=DN(RN)+2 

1C6 C FNF=PN IF TM=0 ANO SUM=0# IF tTM/TMAX t * *2 ( =TMM2 ) 

107 C =0 OR IF !SUM/TMAX!**2(=SUMM2) .NE. AND 

108 C FEM2(=TVM2/SUMM2) .LT. T0LER2 
1C9 C CONTINUED FF-ACTICN 

110 C AZR .CT. 1. RHOLS .LE. RHO .LT. RHOLC 

HI c RHOLC=S0RT(ONE/(TWC*TOLER) ) 

112 C EXF(AZ**2)*ERFC(AZ)= 

113 c <2*AZ/SQRT(PI >)*< 1 I/I { 2*( AZ**2 ) + l )- 

114 C 1*2 I/I <2*<AZ**2l+5>- 

115 C 3*4 I/I (2*(AZ**2)+9)-...) 
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= <£*AZ/PTPI) *I I( AM(RN) I/I EV(RN)) 

PN=1.2.-...RNF 
A*< 1 )=1 

*f(tJNH J=-*M<RN+1 )*<WMRN+1) + 1) 
P*( 1>=2*<AZ**2»+1 
3M(RK*1 )=E«(RN)+4 
WM{ I )=-l 

*M( f-N+1 J=WM(RN) *2 
=( AZ*(FM/GM) )*2/RTFI 
=<AZ»F <PN>)*2/PTPI 
F W { - 1 ) = 1 
GK(-l )=0 

Ff<aj=a 

GV(0)=1 
FC (KNJ=BM{HNJ#FM(RN-1 )+AM<RN) *FK(Rn'-2 J 
GM«N) = BM<RN J*GM(FN-l ) + AM (RN ) * GM SN-2 1 
FNF=PN IF KEM2(F0R J*. E .= ( F t RN )-F ( RN- 1 ) 1 /F (RN) ) 

• LT. TCL2 OR 
KNF = frN-l IF R£M2(RN> . GE. REM2(PN-1) 
ASYMFTCTIC EXFANSICN 

AZF< .LE. 1. AEl_L .LE. RFC .LT* RHOLC 
(FCR FKLC .LE. RHC .LE. CVAX. TO PRESERVE 
ACCURACY" AN ALTERNATIVE CGMPUTAT ICN QF 
THE FIRST TERM CF THE A.E. IS EMPLOYED.) 
EXF(AZ**£)#eFFC(AZ)=(5UM(5GN{FN)*TM(RN)> )/SQRT(PI I 

RN=Q, 1, .. <( RNF 
SGN(0)=1 

SGN(FN+l)=-SGN(RN) 
TK(PN )=( 1/AZ)*{ 1»3.. . (2*RN-1 ) )/ 

(2*< AZ**2) )**RN 
TM(C1=1/AZ 

TMRN*1)=CN(KN+1 )*( TM ( RN )*(!/( 2* (AZ**2)) i) 
CM 1) = 1 

CMFN + l)=DN(RN)+2 
FNF=PN IF PEM2(FGR R . E.-TK/SUM ) .LT. T3L2 
PNF=FN-1 IF TWM2CRN) .GE. TMPK2 ( = T.MM2 (RN- 1 J ) 

(DIVERGENCE) 
RANGE. 

fcXP(Z**Z)i>EPFC{2) IS VALID FOR ZR .GE. THRCUGHOUT 
THE EMI« "ACMNE RANGE, ERF ( 2 ) . ERFC (Z > AND 
EXF(Z**2)»ERFC(Z)(FCR ZR .LT. C) ARE LIMITED EY THE 
RANGE *N0 ACCURACY OF TFE SINE. COSINE AND/OR THE 
EXFCNeNTlAL LIERASY ROUTINES. 
ACCURACY. Tl-c MAXIMUM RELATIVE ERRCC (GENERALLY IN ERFC 1 
-XCEPT IN TFE IMMEDIATE NEIGHECRHCCD CF ZERCSt 

is e<-e) cn ti-e umvac noe FCR 

SINGLE FPECISICN CLMPUT AT ICN. THE REAL 
AND IMAGINARY PARTS INDEPENDENTLY AS *ELL 
A£ TFEIF ZEPCS ENTER INTC CONS IDEP ATI CN . 
PRECIS1CN. VARIABLE - EY SETTING A PREDETERMINED VALUE OF 

TCLER 
MAXIMUM LNIVAC 1106 TIME/SHARING EXECUTIVE SYSTEM 
TIMING. S.F. D.P. 

(SECCNCS) .0101 .052 

STCRAGE. 1171 KCRCS REQUIRED EY THE UNIVAC 110E CCMPILER 
(212 FCRTRAN STATEMENTS. 95 VARIABLES} 
C* TFE PFCRT VER I F I ER , A . C.HALL AND B.G.RYDER 

(HELL LAECPATCRIES, MURRAY HILL. N.J.J PRCC-CF THE COMPUTER 
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174 C SCIENCE AND STATISTICS EIGHTH ANNUAL SYMPOSIUM ON THE 

175 C INTERFACE, UMV. CF CALIF. . LOS ANGELES* FEE. 13-1 A. 1975. 

176 c***** 

177 SUBROUTINE ERRZ< 2R.2I . ERFZR .ERFZI .ERFCZR.ERFCZI t 

178 1 EZ2CZR.EZ2CZI, IERR) 
175 C MACHINE CEFENDENT CONSTANTS 

180 CMAX=. 1701411 E2E39 

181 CMIN=. 14693C754E-38 
ie2 ULSC=. 104857eE7 

183 T0LER=.745E-8 

164 C NCTE TOLER IS SET TC THE PRECISION CF 

165 C THE UNIVAC 1108 SINGLE PKECISIGN ARITHMETIC. 

166 RTPI=1.77245365090E516027DO 

167 C CThER CCr»STANTS 
iee ZERC=0 

165 CNE=1 

ISO T*C=2 

191 THPEE=2 

152 F0UF=4 

153 CNPTFV=THREE/T*C 

154 C INITIALIZATION OF ERROR INDICATORS 

155 IERR=0 

156 IQ=0 

157 C FLNCTICN REFERENCES 

158 C NOTE FUNCTICN REFERENCES OCCUR IN THE REGIONS OF STATEMENT 
155 C LABELS 5.15 AND 85 AND IN STATEMENT LABELS U0.31S.S1S. 
20C C 517 AND 95£. 

201 C SET UP FCR Z IN FIRST QUADRANT AZ = AZR+ I AZ I 

202 5 AZR=AES(ZR) 
20J AZI=AES(ZI) 

2C4 ARIfN=AMINl < AZF .AZI ) 

205 AKIVXsAMAXl <A2R ,AZI) 

2CC IF (ARIMX .GT. ZERO) GO TO 10 

2C7 C SFECIAL CASE. Z=0 

2C6 ERFZR=2ER0 

20^ EPFZI=ZEFO 

210 EPFCZP=CNE 

211 ERFCZI*Z£RC 

212 EZ2CZF=CNE 
2\Z EZ2CZI=ZERC 
2 14 RETURN 

215 C CCNTROL VARIAELES 

21f 10 T0LER2=TCLEfc*TCLER 

£17 TCL2=TCLEP2/FCLR/TWC 

21c r?HCLS = CNPTFV 

215 RMM*X=AFIVN/AFIMX 

220 15 CMAXLN=ALCG(CV4X) 

221 CMINLN=ALOG(CMN) 

22? AELL=SCRT(-ALCG<TCLER) ) 

223 3HCLC = SCRT(CNE/(TWC*TOLER ) ) 

224 PRHC=SGRT ( RMN NX*KMNMX+CNE ) 

225 C CCMPUTATUN OF AUXILIARY QUANTITIES 

226 C CCMFLTATICN OF AZ**2=Z2R+I AZ2I 

227 IF (AFI*X .LT. CNE) GC TO 60 

226 C CVERFLOW CHECK CN Z2PN=-Z2R 

2 29 TEKP3={ (CNE-RVr MX ) * AR I MX > * (CNE+RMNMX ) 

333 IF (TE"PE .LT. CMAX/ARIMX) GO TC 20 

22 1 TEKPC=CMAX 
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GO TO 20 
20 TEMFC=TEMPE*APIWX 
30 IF (AZP-ARIMX) £0#40.50 
40 Z2RN=-TEHPC 

GC TC 70 
£0 Z2PN=TEMPC 

GO TO 70 
60 Z2RN=(AZI+AZR )*<AZI-AZR J 
GC TC 80 

CVERFLOW CHECK CN AZ2I 
70 IF (ARIMN ,LT. < CMAX/TWG > /ARI MX > GC TO 80 
AZ21=CKAX 
GO TC ?0 
80 AZ21=(TW0*ARIVN>*ARIMX 

PRELIMINARY COMPUTATIONS FOR EXP(AZ»*2) AND 
EXPt-AZ**2) 

CHECK IF VALID ARGUMENT FOR SIN/COS 
IF (AZ2I .GE« UL£C) GO TO 90 
85 C0AZ2IsCCS(AZ2I ) 
SIA22I=SIN(AZ2I) 

EXTENDING RANGE OF EXP ROUTINE 
90 TEKF=Z2PN/THREE 
Z2R=-Z2PN 

CVERFLOW CHECK 
IF (TEMP .LT. CVAXLN) GO TO 100 
EMZ2D3=CMAX 
EZ2R=ZERC 
GC TO 190 

UNDERFLOW CHECK 
100 IF CTEfP .C-T. CKINLN) GO TO HO 
EMZ2C2=ZERC 

EXP(AZ**2) OVERFLOWS (IQ=1) 
IQ=1 

GC TC 160 
110 EWZ2D6=EXP(TEMP/T*0) 
EfZ2D2=EKZ206*EyZ2Dfc 
EZ2C6=CNE/EMZ2C6 
IF <EZ20C .LE. CNEI GO TO 180 
J = l 

PEXP=ZEFG 
TEMF=EZ2D6 
120 IF (TEMP ,C-E. CMAX/EZ2D6) GO TO 150 
TEMF=TEMF*EZ2C6 
J=J+1 

IF <J-5) I20tl20*140 
130 PEXP=T£MP 
GC TO 120 
140 EZ2P=TEMP 
GO TO 190 
150 IF (PEXP) 170,160.170 
160 PEXF=CKAX 
170 EZ2R=C*AX 
GO TO 150 
180 PEXP=E22D6**5 

EZ2R=PEXP*E22C6 
TEZ2=(PEXP*TWC)*EZ2DC 

CCMFLTATION OF RHO 
CVERFLOW CHECK 
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290 190 IF IARIKX .LT. CMAX/PRHOl GO TO 200 

291 RHC=CVAX 
252 GC TC 210 

293 200 RHC*ARIfX*PFHC 

294 C KEThCD SELECTION 

295 210 IF {RHC .LT. RHCLC ) GO TO 220 

296 C IMFRCVE ACCURACY FOR LARGE RHO 

297 FA=AZF/AfiIMX 
29E FB=AZI/ARIMX 

299 FC=RTPI*(RMh^X*RMNKX+ONE) 

300 EZ2CZR*=<FA/ARIMX)/FC 

301 EZ2CZI=-(FB/APIMX)/FC 

302 GO TO eOO 

303 220 IF (RHC .LT. RHCLSJ GO TO 300 

304 IF (AZR .GT. CNE1 GC TO 500 
3C5 IF (RHC-AELL) 300i70Q,700 

306 C PCWEP SERIES FOR ERFiAZJ 

3C7 C INITIALIZATION 

306 300 5U*R=Z£RG 

309 SUMI=ZGFG 

210 SGN=CNE 

311 RN=ZERC 

2 12 ON=CNE 

313 PTMfi=AZR 

314 PTMI=A2I 

21S C COMPUTING SUM 

316 310 TfR=FT*R/DN 

317 TM«FTH/CN 

318 £UMR=SUMR+TVR*SGN 

319 SUKI = SUMI+TWl*SC-N 

320 C SCALING TO AVOIC OVERFLOW OR UNOER- 
221 C FLCW IN APPROXIMATING R.E. 

322 315 TMAX=AMAX1 ( AES (THR ) , AES CTMI ) .ABS(SUMR), ABS(SUMI) ) 

323 IF {TMAX) 320i360i320 

224 320 TVM2=(TKfi/TWAX)**2*<TKI/TMAX)#*2 

225 SUM»'£=<SUMR/TMAX)**2-MSUMI/T!*AX>**2 
326 IF (TVM2) 330i360,230 

227 330 IF (SUVM2) 240.350t340 

323 340 RE*2 = TMM2/SL*<-W2 

229 C TOLERANCE CHECK 

330 IF (REN2 .LT. TCLER2) GC TO 360 

221 C ADDITIONAL TERMS 

232 3S0 Dr*=D^TWC 

323 RN=RK+CNE 

334 SGN=-SGN 

235 TEMF=(FTMR*22fi-PTMI*AZ2l )/RN 

336 PT*»I = f FTMF*AZ2I+PTMI*Z2R >/RN 

337 PTMR=TEMP 

328 GC TC 310 

329 C FUNCTICNS EVALUATED IN FIRST CUADRANT 
340 360 ERFZR=SUKR*TfcC/FTPI 

241 ERFZI=SUMI*TWC/FTPI 

342 EFFCZR=CNE-ERFZP 

343 ERFCZI=-ERFZI 

344 EZ2CZP=(CCAZ2I*ERFCZR-SIAZ2I*ERFCZI >*EZ2R 
3«S EZ2CZI=CSIAZ2I*ERFCZR+CGAZ2I*ERFCZI)*EZ2R 

346 GC TC $40 

347 C CONTINUED FRACTION FOR EXP { AZ**2 ) *ERFC{ AZ J 
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34 1 C INITIALIZATION 

34-.; 500 WV = -CN£ 

3EC SMP=T*C*Z2P+CN£ 

351 BMI=T*C*AZ2I 

352 AM=CNE 
253 FM*2R=CNE 
354 F^2I = ZERC 
36E GMM2R=ZEK0 

356 Gf/M2I = Z£RC 

357 FVM1R=2ER0 

358 FKflI=ZECa 
3E^i GHV1R = CNE 

360 GMVUsZEFO 

361 RERN2=CWAX 

362 FPfi=ZEFC 

363 FPI=ZEPQ 

364 C RECURRENCE RELATION 
3t£ 510 FWR = EfR*F^MlR-BMI*FMMl I-f AM*FMM2R 
366 FWI = EVI*FKVlfi+BMR*FMMl I+AM*FMM2 I 
3 67 GI¥R = ENR*GMM1R-BMI*GMM1 I+AM»GMM2R 
368 GM = EWI*GMWiq + EMR*GMMlI*AM*GMM2I 

365 C CCNVERGENT F=FM/GM 

37C C SCALING TO AVOID OVERFLOW IN 

371 C CCMPUTING CCNVERGENT 

372 £15 TVAX=AMAX1 ( AES ( FMR ) , A6S ( FMI )»ABS(GMR) ,ABS(GMI >> 

373 $FKR=F*R/TMAX 

374 SFMI=FMI/TMAX 
37= SGMR=GMR/TMAX 

376 SGVI = GVI/T/«AX 

377 TEMF = SGfR*SC-MR + SGMI*SGMl 

376 FR=(SFMR*SGH«P4SFMI*SGMI)/TEMP 

3 75 FI=(SFMI*SGVP-SFMR*SGMI)/TEMF 

360 C APPROXIMATING R.E. 

aei TEMF=Ffi*FR+FI*FI 

362 TEMFA=FR~FPR 

383 TEMPE=FI~FFI 

364 REM2 = ( TEMFA*T EMPA+TEMPB* TEMP B ) /TEMP 

365 C TOLERANCE CI-ECK 

366 IF (REM2 .LT. TCL2) GO TO 530 

367 IF (REM2 .GE. REPM2) GO TO 520 

3E8 C ADDITIONAL CONVEPGENTS 

365 H*=*V+T*C 

350 EMR=EMR+FCUR 

351 A*=-*W*( W*+CNE) 

352 FVM2F=FMM1R 

353 FfV2I=FMMlI 

354 GMM2P=G*>M1R 

355 C-NV2I = GKV1I 

356 FMM1R=FMP 

357 FfVlIsFKI 
35S GMM1R=GMR 
355 GYV1I=GMI 
4C0 fpr=fr 

401 FPI=FI 

402 REFf£=REM2 

403 C SCALING 

404 C SCALING SHCULD NCT EE DELETED AS TFS VALUES OF FMH.FMI 

405 C CMR.GMI WAV CVERFLC* FCR SMALL VALUES OF REAL OF Z 
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40* £17 AB« = AVAX1 <AE c (e.MP KAdS(PMI) ) 

4S7 if (TV4X .LT. (CMAX/FCUK>/(TV(C*6BM-AM) ) GC TO 510 

40t f^^Ssfn^z^/tkax 

4C FMM2 I=FfM2I/TWAX 

410 G*M2PsG*'W£P/TVAX 

411 0t-V2 I=GWK2 I/TMAX 

412 FVtf lp=Ff>«lR/TMAX 
4 1 3 F^MI I = FV^1 I/TWAX 
41<i GVV 1P=GN'«1K/TMAX 
MS- GWV l I = Gf'f/l I/TWAX 
MS GC TC e 10 

417 C RELATIVE ERROR I NCRE ASED-RQUNDCFFS 

41* C ACCEPT PRIOR CCNVERGENT 

dl'j t20 FS=FPP 

«20 FI=FPI 

421 C EVALUATE EXP ( AZ **2 ) *ERFC ( AZ 1 

422 120 EZ2CZF=IAZR*FK-AZ1*F1)*TKC/RTPI 
42 J EZ2CZI=(AZI*FF+AZR*FI)*T*C/RTPI 

424 GC TC SCO 

425 C ASYVPTCTIC EXPANSION FCR EXP( AZ** 2 ) *ERFC ( AZ ) 
42t C INITIALIZATION 

427 700 TZ2R=T»C*Z2* 

428 TZ2I=TkC*AZ2I 

42^ TEVF=TZ2R*TZ2F+TZ2I*TZ2I 

420 RTZ2R=TZ2R/TEyP 

421 RTZ2I=-TZ2I/TEWF 

432 TMM1R=(AZR/RHC)/RHC 

433 TMMll=-(AZI/fcHC)/RFC 

4 34 TMFV2=TMMlR*T^y Ifi+T^Ml UTMMJI 

4cE SUVRsTI^WlR 

42£ SUMI=TWflI 

4?7 ONsCNE 

43E SGr=-Cf>E 

42<3 C CCMFUTING SUM 

440 710 TVR = CN*(T»'MIR*RTZ2R-TMMI I*RTZ2I) 

44 1 TVI=CN*(TMM1I*RTZ2R+TMM1R*RTZ2I ) 

442 SUMR=£GN*TMP+SUMR 

443 SUMI = SGN*TM + SUMI 

444 C APPROXIMATING R.E. 

445 SU^2sSUMR*SL*R + SUMl*SUMI 
44C TCf'2 = TWK*T^F + TWI*TMI 

447 REN2=TPM2/StVM2 

44e C TCLERANCE CFECK 

449 IF <REM2 .LT. TCL2> GO TO 730 

4EC IF (TVM2 .LT. TMPK2) CO TO 720 

4E1 C DIVERGENT PATH 

4*2 SUMBsSUMR-SGN*TMR+SGN*T^MlR 

4*3 SLMI=SUMI-SGK*TMI*SGN*TMM1I 

4£4 GC TC 730 

455 C ADDITIONAL TERMS 

4S<5 720 SGN=-SGN 

4S7 Or* = DtH-Tl»0 

4Ee TVM1F=TMR 

459 T^MlIsTMI 

460 TMFW2 = Tf'M2 

461 GC TC 710 

4*2 C EVALUATE EXP( AZ**2 ) *ERFC ( AZ ) 

4C3 730 EZ2CZR=SUMR/RTPI 
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4<4 EZ2CZ I = SUMI/RTPI 

4*5 C MAINTAINING ACCURACY IN EZ2CZR 

46e C FCR SMALL AZR 

4C7 TEMF=E22CZR*TCLER 

4ee IF (TEMP .GT. EZ2R) GO TO 800 

<i( c . TEMF = CHN*FHCLC 

470 IF (EZ2CZR .GT. TEMP) GO TO 750 

471 IF (EZ2R .GT. TEMP) GO TO 750 
47? IF (AZF) 74C.7EC.740 

473 C INDICATE RESULTING ERRORS IN ERFC 

474 C (AND ERF)(IERR=2) 
47? 740 IERfi=2 

476 750 EZ2CZR = EZ2CZfc-f EZ2F*C0AZ2I 

477 C EVALUATE ERFC AND ERF FOR AZ IN CONTINUED 
47E C FRACTION AND ASYMPTOTIC EXPANSION REGIONS 

475 800 IF (A22I .LT. ULSC) GO TO 830 

4£0 C INVALID ARGUMENT FOR SIN/COS 

4°1 C ERFC(AND ERF) INVALID (IERR=2) 

462 IF (AZF-ARIMX) 620.eiO.820 

432 C AZR .GE. AZI 

4E4 810 ERFCZR=ZERC 

4£ C J ERFCZI=ZERC 

nee GC TC S10 

4F7 C AZR .LT. AZI 

4ee e20 erfczr=cvax 

489 ERFCZI=CKAX 

450 GC TC S10 

451 C VALID ARGUMENT FOR SIN/COS 
4?r 830 IF (RFC .GE. RHCLC) GO TO e40 

4^3 C RHC .LT. RHOLC 

4S4 TEMPA=EZ2CZR*CCAZ2I+EZ2CZ I*SIAZ2I 

4CS TEMFB=-EZ2CZR*SIAZ2I+EZ2CZI*COAZ2I 

4St TEMPC=EWZ2D3 

4<?7 GC TC E50 

4C t c Rt-C .GE. RHOLC 

4S<; 640 TFMFA=(FA*CCAZ2I-FE*SIAZ2I )/FC 

5C0 TEMPB=<-FA*SIAZ2I-FB*COAZ2I )/FC 

SCI TEVF=ENZ2D3/AFiyX 

50< TEMFA=TEMPA*TEKF 

5.03 TEMFE=TEMPE*TEVF 

504 TFMFC=CNE 

5C3 F60 IF (EVZ233 .LE. CNE ) GO TO 920 

50f C EVALUATE ERFC(AZ)(AZI . GT . AZR) 

507 1=1 

5CF TEMF=TE^FA 

50«, eeo J = l 

eio SGN=CN5 

Ell FC=TEMFC 

Oi£ IF (TEMP .CE. ZERO GO TO 870 

513 SGN=-SGN 

fl4 TEMF=-TEMP 

fl^ g70 IF (TEMP .LT. CMAX/FD) GO TO 880 

hit TENF=CMAX 

517 IC=2 

He GC TG e?0 

j.19 fc«0 TE^R=TEMF*FO 

5.3 -j J = J+1 

t £! IF (J .C-T. 3) GC TO 890 



679 



522 FD=EMZ2C3 

£23 GC TO 67C 

£24 8QO IF (I .GE. 2J GO TO 900 

£25 ERFC2R=TEXF*SCN 

£26 1=1+1 

£27 TEMF=TEMPE 

£28 GC TC 860 

£29 900 ERFCZI=TEMP*SGN 

£30 IF CIQ .NE. 2) GO TO 930 

£31 C ERFC(AND ERF) INVALID C IERR=2 ) 

532 C OVERFLOW CF ERFC ( AZ > ( 10=2 > 

£32 910 IEPR=2 

£24 GC TC «30 

£35 C EVALUATE ERFC(AZ)(AZI .LE. AZR) 

E3€ 920 ERFCZR=(TEMPA*EMZ202)*EMZ2D3*TEMPC 

£37 ERFCZI=(TEMPE*EKZ2D3»*EMZ2D3*TEMPC 

S28 C SPECIAL CASE (AZR=0) 

£29 930 IF (AZR .LE. 2EPG) ERFCZR = CNE 

E4C C EVALUATE ERF(AZ) 

£41 ERFZR=CNE-ERFCZR 

£42 ERFZI=-ERFCZI 

543 C SYMMETRY RELATIONS APPLIED 

£44 940 IF (ZP-AZR1 950,1000.950 ' 

£45 C REAL CF Z .LT. 

E4C 950 ERFZR=-ERFZR 

£47 ERFCZR=T*0-ERFCZR 

£48 IF (A22I .GE. ULSC) GO TO 960 

£49 IF (IC #EQ. 1) GO TC 960 

550 C KAIf>TAINlNG ACCURACY IN 2*EXP(Z**2) 

£51 IF (A2R .LE. #21) GO TO 980 

££2 9EE TEMP = AI»AX1 ( AES { S I AZ2 I ) . ABS(C0AZ2in 

£53 IF (TEMP .LT. ( ( CM AX/TW0)/PEXP } /EZ2D6) GO TC 970 

£54 C EXP(Z**2)*ERFC(ZI INVALID (IEPR+1) 

£55 96C IERR=IERR+1 

556 EZ2CZP*CMAX 

£E7 EZ2CZI=CVAX 

£5£ GO TO 1000 

5E9 970 TEZ2R=((PEXP*CCAZ2I>*EZ2D6J*T*0 

5C0 TE22Is((PEXP*SIAZ2I)*EZ2D6)*TWC 

5«1 GC TO 990 

£€2 980 TE22R=TEZ2*CCAZ2I 

£C3 TEZ2I=TEZ2*SIAZ2I 

5*4 C EVALUATE EXPC Z** 2 ) *ERFC ( Z ) 

£C£ 990 EZ2CZRSTEZ2F-EZ2CZR 

5C6 EZ2CZI*EZ2CZI-TE22I 

St7 1000 IF (Zl-AZI) 1010.1020.1010 

5C8 C IMAGINARY CF Z .LT. 

ttc 1010 ERFZI=-ERFZI 

57C ERFCZU-ERFCZI 

£71 EZ2CZI=-EZ2CZI 

£72 1020 RETURN 

£72 END 
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TABLE 1 
























ERFC(Z) 














RHO\THETA 




15° 


30° 




37.5° 


45° 




.00 


.1000004-01 


.0 


.100000+01 


.000000 


.100000+0 1 . 


000000 


.100000+01 . 


000000 


.100000+01 


.000000 




.02 


.977435*00 


.0 


.978204+00 ■ 


-.583879-02 


.980456+00 -. 


112808-01 


.982095+00 -. 


137355-01 


.984040+00 


-.159556-01 




.04 


.954889+00 


.0 


.956420+00 • 


-. 116648-01 


.960912+00 -« 


225435-01 


.964183+00 -. 


274543-01 


.968068+00 


-.318984-01 




.06 


.932378+00 


.0 


.934662+00 • 


-.174654-01 


.941368+00 -. 


337702-01 


.946257+00 -. 


411398-01 


.952070+00 


-.478156-01 




.08 


.909922+00 


.0 


.912942+00 • 


-.232279-01 


.921824+00 -. 


449428-01 


.928310+00 -. 


547751-01 


.936033+00 


-.636943-01 




.10 


.887537+00 


.0 


.891273+00 • 


-.289397-01 


.902280+00 -. 


560434-01 


.910337+00 -. 


683437-0 1 


.919946+00 


-.795217-01 




.12 


.865242+00 


.0 


.869667+00 • 


-.345886-01 


.882738+00 -. 


670542-01 


.892330+00 -. 


818288-01 


.903796+00 


-.952846-01 




.14 


.843053+00 


.0 


.848138+00 - 


-.401625-01 


.863197+00 -. 


779575-01 


.874283+00 -. 


952136-01 


.887571+00 


-.110970+00 




. 16 


.820988+00 


.0 


.826697+00 ■ 


-.456493-01 


.843657+00 -. 


887357-01 


•856190+00 -. 


108481+00 


.871258+00 


-. 126564 + 00 




.18 


.799064+00 


.0 


.805358+00 


-.510376-01 


.824121+00 -. 


993713-01 


.838045+00 -• 


121615+00 


.854845+00 


-. 142053+00 




.20 


.777297+00 


.0 


.784132+00 ■ 


-.563161-01 


.804590+00 -< 


109847+00 


.819844+00 -. 


134598+00 


.838321+00 


-.157424+00 




.30 


.671373+00 


.0 


.680129+00 


-.806921-01 


.707071+00 -« 


159241+00 


.727824+00 -. 


196662+00 


.753652+00 


-.231995+00 




.40 


.571608+00 


.0 


.580740+00 


-.100871+00 


.610081+00 -. 


202202+00 


.633845+00 -. 


2S2418+00 


.664672+00 


-.301347+00 




.50 


.479500+00 


.0 


•487320+00 


-.1 15988+00 


.514270+00 - 


.236937+00 


.537866+00 -. 


299764+00 


.570447+00 


-.363359+00 




.60 


.396144+00 


.0 


.401058+00 


-.125565+00 


.420627+00 -. 


261983+00 


•440342+00 -. 


336667+00 


.470518+00 


-.415645+00 




.70 


.322199+00 


.0 


.322926+00 


-.129533+00 


.330478+00 - 


.276308+00 


.342317+00 -i 


361274+00 


.365070+00 


-.455585+00 




.80 


.257899+00 


.0 


.253631+00 


-.128213+00 


.245436+00 - 


.279418+00 


.245465+00 -. 


372052+00 


.255106+00 


-.480414+00 




.90 


.203092+00 


.0 


.193577+00 


-.122258+00 


.167296+00 -. 


271437+00 


.152078+00 -i 


367979+00 


.142621+00 


-.487404+00 




1 .00 


.157299+00 


.0 


. 142847+00 


-.1 12577 + 00 


.978859-01 - 


.253165+00 


.649787-01 -« 


348749+00 


.307358-01 


-.474148+00 




1 .10 


.119795+00 


.0 


.101201+00 


-.100228+00 


.388720-01 - 


.226077+00 


-.126659-01 -, 


314972+00 


-.762502-01 


-.438939+00 




1 .20 


.896860-01 


.0 


.680972-01 


-.8631 11-01 


-.844458-02 - 


. 192259 + 00 


-.776217-01 -, 


268341+00 


-.172938+00 


-.381252+00 




1 .30 


.659921-01 


.0 


.427427-01 


-.718743-01 


-.433244-01 - 


.154271+00 


-.126958+00 -• 


211694+00 


-.253111+00 


-.302268+00 


o\ 


1 .40 


.477149-01 


• 


•241557-01 


-.578252-01 


-•657413-01 - 


•114934+00 


-•158478+00 -« 


148963+00 


-•310259+00 


-•205369+00 


00 


1 .50 


.338949-01 


• 


.112438-01 


-.448760-01 


-.764609-01 - 


.770708-01 


-.171150+00 -i 


849352-01 


-.338390+00 


-.965018-01 


H-* 


1 .60 


.236516-01 


• 


•288326-02 


-.335145-01 


-.770202-0 1 - 


.432136-01 


-.165475+00 -, 


248347-01 


-.333094+00 


. 15760 0-01 




1 .70 


.162095-01 


.0 


-.200643-02 


-.240031-01 


-.696000-01 - 


•153368-01 


-.143694+00 . 


262601-01 


-.292775+00 


. 120601+00 




1 .80 


.109095-01 


.0 


-.44 0212-02 


-.164025-01 


-.567976-01 


.535338-02 


-.109748+00 


.641075-01 


-.219843+00 


.205946+00 




1 .90 


.720957-02 


.0 


-.513475-02 


-. 106123-01 


-.413301-01 


.185261-01 


-.689521-01 


.859866-01 


-. 121580+00 


.260 096+00 




2.00 


.467773-02 


.0 


-.486777-02 


-.641967-02 


-.257133-01 


.247077-01 


-.273484-01 


.912368-01 


-.103117-01 


.273926+00 




2.10 


.297947-02 


.0 


-.409504-02 


-.354979-02 


-.119733-01 


.251347-01 


.916428-02 


815039-01 


.974865-01 


.243355+00 




2.20 


.186285-02 


• 


-.315468-02 


-.170957-02 


-. 144 281-02 


.215037-01 


.358034-01 


.605788-01 


.183390+00 


. 171536 + 00 




2.30 


.114318-02 


• 


-.225377-02 


-.621965-03 


.531949-02 


.156622-01 


.497620-01 


.337876-01 


.230486+00 


.700421-01 




2.40 


.688514-03 


• 


-.149843-02 


-.481389-04 


.847425-02 


.930477-02 


.507367-01 


.700451-02 


.227631+00 


-.417236-01 




2 .50 


.406952-03 


.0 


-.924217-03 


.201683-03 


.872881-02 


.372996-02 


.409438-01 - 


•145264-01 


.173531+00 


-. 139462 + 00 




2.60 


.236034-03 


.0 


-.523179-03 


.266848-03 


.709056-02 - 


.293708-03 


.245557-01 - 


•272819-01 


.792740-01 


-. 199249+00 




2.70 


.134333-03 


.0 


-.265334-03 


.241454-03 


.460569-02 - 


.253794-02 


.664754-02 - 


•301636-01 


-.321336-01 


-.204295+00 




2.80 


.750132-04 


.0 


-.113783-03 


.183127-03 


.214179-02 - 


.321372-02 


-.809354-02 - 


•245881-01 


-.130194+00 


-. 151392+00 




2.90 


.410979-04 


.0 


-.339435-04 


.122861-03 


.256793-03 - 


.279917-02 


-.165397-01 - 


.138997-01 


-.185124+00 


-.545715-01 




3.00 


.220905-04 


.0 


.198248-05 


.741051-04 


-.830403-03 - 


.184961-02 


-.178141-01 - 


•225516-02 


-. 178018+00 


.564096-01 




3.50 


.743098-06 


.0 


.383634-05 


-.331302-06 


.361438-04 


•343118-03 


•663124-02 


.774278-03 


.146013+00 


-.667835-01 




4.00 


.154173-07 


.0 


-.506480-07 


-.121829-06 


-.10C583-04 - 


.454411-04 


-.206887-02 


.808579-03 


-.704373-01 


. 121816+00 




4.50 


.196616-09 


.0 


-.172664-08 


.241877-08 


•342245-05 


.358814-05 


.149447-03 - 


.641621-03 


-.697054-01 


-. 103987+00 




5.00 


.153746-11 


.0 


.431344-10 


-.812606-11 


-.410390-06 


.690194-07 


.163226-03 


.593711-04 


.909031-01 


-.666628-01 




5.50 


.735785-14 


.0 


-.401550-12 


-. 137042-12 


-.988939-10 - 


.274562-07 


.529439-06 


.406179-04 


.945377-01 


.396363-01 




6.00 


.215197-16 


.0 


.222076-14 


. 150084-14 


.136924-08 - 


•383336-09 


-.550781-05 


•636094-05 


.563713-01 


.752047-01 




6.50 


.384215-19 


.0 


-.905982-17 


-.633184-17 


.477805-10 


.324007-10 


-.126107-05 


.886874-06 


.498547-01 


.710156-01 




7.00 


•418383-22 


.0 


.276523-19 


. 109112-19 


•941894-12 


.157585-11 


-.164519-06 


.187753-06 


.71085e-01 


.379415-01 




7.50 


.277665-25 


.0 


-.517996-22 


.562585-23 


.233092-13 


.393021-13 


.525818-10 


.356932-07 


.667858-01 


-.345871-01 




8.00 


.112243-28 


.0 


.396093-25 


-.443759-25 


.731135-15 


.506747-15 


.420606-08 


.160997-08 


-.258232-01 


-.656143-01 




8.50 


.276232-32 


• 


• 11 1560-28 


.427748-28 


•130972-16 - 


.343079-17 


.123912-09 - 


.485650-09 


-.475602-01 


.462887-01 




9.00 


.413703-36 


.0 


-•213063-31 


-.174915-32 


.313042-20 - 


. 160996-18 


-.462157-10 


.168052-10 


•623821-01 


-.612171-02 



RHQVTHETA 

.00 .100000 

.02 .965491 

.04 .970967 

.06 .95641 1 

.06 .941809 

.10 .927144 

.12 .912401 

.14 .897565 

.16 .862621 

.18 .867552 



50 
♦01 .000000 
♦00 -.172863-01 
♦00 -.345635-01 
♦00 -.518226-01 
♦00 -.690545-01 
♦00 -.862497-01 
♦00 -.103399+00 
♦00 -.120493+00 
♦00 -. 137521+00 
♦00 -.154473+00 





.20 


.852345+00 


-. 171339 + 00 




.30 


.773712+00 


-.253982+00 




.40 


.669468+00 


-.332643+00 




.50 


.59e041 +00 


-.405420+00 




.60 


.498195+00 


-.469950+00 




.70 


.389205+00 


-.523342+00 




.80 


.271098+00 


-.562163+00 




.90 


. 144933+00 


-.582512+00 


1 


.00 


.131385-01 


-.5e021 1 +00 


1 


.10 


-.120165+00 


-.551 166+00 


1 


.20 


-.248682+00 


-.491928+00 


1 


.30 


-.364829+00 


-.400481+00 


1 


.40 


-.457890+00 


-.277246+00 


On 1 


.50 


-.516623+00 


-.126215+00 


00 , 


.60 


-.529455+00 


.439353-01 


.70 


-.486589+00 


.219059+00 


i 


.80 


-.382619+00 


.379466+00 


i 


.90 


-.219614+00 


.50 1 C81+00 


2 


.00 


-. 101646-01 


.558267+00 


2 


.10 


.220530+00 


.S2e666+00 


2 


.20 


.434783+00 


.399806+00 


2 


.30 


.5861 12+00 


.176582+00 


2 


.40 


.627560+00 


-.112358+00 


2 


.50 


.524677+00 


-.412231+00 


2 


.60 


.271303+00 


-.647141+00 


2 


.70 


-.961516-01 


-.735157+00 


2 


.80 


-.492661+00 


-.61 3686+00 


2 


.90 


-.794447+00 


-.271241+00 


3 


.00 


-.e701 18+00 


.227570+00 


3 


.50 


.128375+01 


-.438817+00 


4 


.00 


-. 143233+01 


.176970+01 


4 


.50 


-.153046+01 


-.394580+01 


5 


.00 


.818979+01 


-.289861 +01 


5 


.50 


. 141 149+02 


.136673+02 


6.00 


.892289+01 


.480481+02 


6.50 


.767743+01 


.133290+03 


7 


.00 


. 164446 + 03 


.364872+03 


7, 


50 


.126058+04 


.375323+03 


a 


00 


.231 199+04 


-.413300+04 


a 


.50 


-.181429+05 


-.441024+04 


9 


.00 


.404091 +05 


.697050+05 



.100000+01 
.988713+00 
.977408+00 
.966067+00 
.954672+00 
.943204+00 
.931646+00 
.919978+00 
.908183+00 
.896242+00 

.884135+00 

.820454+00 

.749698+00 

.669242+00 

.576286+00 

.467898+00 

.341 135+00 

.193283+00 

.222736-01 

-.172655+00 

-.389989+00 

-.624572+00 

-.865837+00 

-.109565+01 

-.1 28600+01 

-.139708+01 

-.137659+01 

-.1 16196+01 

-.687539+00 

.100155+00 

. 121649+01 

.260024+01 

.406480+01 

.524907+01 

.559076+01 

.436233+01 

.827238+00 

-.542122+01 

-. 139065+02 
.436770+02 

-.2e6672+03 
.303475+04 

-.236016+05 

-.192298+06 
.434927+07 
.13011 1+09 
.3C8919+10 
. 107926+12 
.557097+13 
.231401+15 

-.1 17844+17 



60 

.000000 
-. 195441-01 
-.390882-01 
-.586322-01 
-.781761-01 
-.977195-01 
-. 1 17262 + 00 
-.136803+00 
-. 156342+00 
-.175878+00 

-. 155409+00 

-.292919+00 

-.389843+00 

-.485367+00 

-.578076+00 

-.665719+00 

-.744933+00 

-.810912+00 

-.857061 +00 

-.874670+00 

-.852699+00 

-.777828+00 

-.635009+00 

-.408856+00 

-.863300-01 

.338770+00 

.85901 1 +00 

. 144402 + 01 

.202914+01 

.250402+01 

,270506+Cl 

.241951+01 

.141207+01 

-.51 1462+00 

-.339050+01 

-.694257+01 

-. 1 03663+02 

-. 121876+02 

-.102944+02 

.609546+C2 

-.315634+03 

.899044+03 

.194203+05 

-.331209+06 

-.444020+07 

.950077+07 

. 172358+10 

.606848+1 1 

.462008+12 

-.228725+15 

-.213680+17 



ERFC(Z) 








70 




80° 


.100000+01 


.000000 


.100000+01 


.000000 


.992279+00 - 


.212081-01 


•996080+00 


-.222273-01 


.984542+00 - 


.424252-01 


.992150+00 


-.444703-01 


.976774+00 - 


.636604-01 


.988203+00 


-.667446-01 


.968559+00 - 


.849226-01 


.984228+00 


-.890659-01 


.961060+00 - 


.106221 +00 


.980217+00 


-.111450+00 


.953123+00 - 


. 127564+00 


.976160+00 


-.133913+00 


.945070+00 - 


. 148961+00 


.972048+00 


-. 156471+00 


.936906+00 - 


. 170421+00 


.967870+00 


-. 179140+00 


.928612+00 - 


. 191952 + 00 


.963617+00 


-.201936+00 


.920173+00 - 


.213564+00 


.959279+00 


-.224877+00 


.875153+00 - 


.323124+00 


.935924+00 


-.342344+00 


.823615+00 - 


.435932+00 


.908660+00 


-.466099+00 


.762712+00 - 


.552889+00 


.875613+00 


-.598673+00 


.688562+00 - 


.674665+00 


.834330+00 


-.742997+00 


.598036+00 - 


.801524+00 


.781501+00 


-.902533+00 


.484501+00 - 


.933057+00 


.712565+00 


-. 108143 + 01 


.341532+00 - 


.106775+01 


.621 138+00 


-.128468+01 


.160625+00 - 


.120233+01 


.498184+00 


-.151831+01 


-.686592-01 - 


. 133079 + 01 


.330784+00 


-. 178947+01 


-.358812+00 - 


. 144298+01 


.100318+00 


-.210652+01 


-.724161+00 - 


. 152248+01 


-.220245+00 


-.247879+01 


-.118000+01 - 


. 154380 + 01 


-.670390+00 


-.291589+01 


-.174051+01 - 


146840+01 


-. 130809+01 


-.342593+01 


-.241452+01 - 


. 123975 + 01 


-.221880+01 


-.401186+01 


-.319777+01 -, 


777391+00 


-.352885+01 


-.466405+01 


-.405945+01 


291660-01 


-.542509+01 


-.534607+01 


.452009+01 


132454+01 


-.818351+01 


-.596797+01 


.561710+01 . 


328529+01 


-.122098+02 


-.633691+01 


-.585370 + 01 


.610278+01 


-. 180944 + 02 


-.606686+01 


-.512891 + 01 


993582+01 


-.266833+02 


-.441546+01 


•265168+01 . 


148074+02 


-.391538+02 


.862691-02 


.274334+01 


204062+02 


-.570679+02 


.976572+01 


.126293+02 


257378+02 


-.8231 62+02 


.294995+02 


.288616+02 . 


285642+02 


-.116747+03 


.675244+02 


.530807+02 . 


245e83+02 


-. 161012 + 03 


i 136532+03 


.854677+02 . 


644245+01 


-.21 1587+03 


.268092+03 


.122167+03 -. 


371981 +02 


-.253771+03 


.499853+03 


.150602+03 -. 


121463+03 


-.246068+03 


.906353+03 


.185185+04 -. 


696951+03 


.126726+05 


. 1 1001 1+05 


.142504+05 . 


268505+05 


.4141 64+06 


-.267047+06 


.581405+05 -. 


694896+06 


-. 104472 + 08 


-.21 1627+08 


.301677+05 . 


237891+08 


-. 159723+10 


.902680+09 


.279467+09 -. 


117057+10 


. 155714+12 


. 170313+12 


.640472+11 . 


633707+11 


.196750+14 


-.425327+14 


.952019+13 


967322+12 


-.152018+17 


.205488+16 


.442269+15 -. 


156615+16 


.619376+19 


.519553+19 


.371118+18 . 


125665+16 


-.143772+22 


-.670033+22 


.785527+19 . 


138783+21 


-.257275+25 


.897512+25 


.624892+23 . 


370071+23 


.1151 66+29 


-. 168785 + 29 


.554262+26 -. 


688648+25 


-.488030+32 


.526592+32 



.1 + 01 


.000000 


.1+01 


-.225706-01 


.1+01 


-.451593-01 


.1 + 01 


-.677841-01 


.1+01 


-.904633-01 


.1+01 


-. 113215+00 


.1 + 01 


-. 136058+00 


.1+01 


-. 15901 1+00 


.1+01 


-. 182093+00 


.1+01 


-.205323+00 


.1+01 


-.228721+00 


.1 + 01 


-.348949+00 


.1+01 


-.476625+00 


.1+01 


-.614952+00 


.1 + 01 


-.767853+00 


.1 + 01 


-.940283+00 


.1+01 


-. 113867 + 01 


.1+01 


-. 137154 + 01 


.1+01 


-. 165043 + 01 


.1 + 01 


-. 199117+01 


.1+01 


-.241591+01 


.1 + 01 


-.295609+01 


.1 + 01 


-.365696+01 


.1 + 01 


-.458473+01 


.1 + 01 


-.583773+01 


.1 + 01 


-.756418+01 


.1+01 


-.9991 12+01 


.1 + 01 


-. 134718+02 


.1+01 


-. 185648+02 


.1+01 


-.261677+02 


.1+01 


-.377471+02 


. 1 + 01 


-.557397+02 


.1+01 


-.842631+02 


. 1 + 01 


-. 130396 + 03 


.1+01 


-.206519+03 


.1+01 


-.334671+03 


.1 + 01 


-.554777+03 


.1+01 


-.940470+03 


.1+01 


-. 162999 + 04 


. 1 + 01 


-.352823+05 


.1+01 


-. 129696 + 07 


.1+01 


-.801975+08 


.1 + 01 


-. 829827+10 


.1+01 


-. 143210+13 


.1+01 


-.41 1275+15 


. 1 + 01 


-. 196225+18 


.1+01 


-. 155349 + 21 


.1+01 


-.203882+24 


. 1+01 


-. 443245+27 


.1 + 01 


-. 159530+31 


.1 + 01 


-.950078+34 



TABLE 2 



EXP(Z**2) *ERFC(Z) 





RHOMHETA 






.00 


. 100000 + 01 


.0 




.02 


.977826+00 


.0 




.04 


.956418+00 


.0 




.06 


.935741+00 


.0 




.08 


.915764+00 


.0 




.10 


.896457+00 


.0 




.12 


.877791+00 


.0 




.14 


.859740+00 


.0 




.16 


.842277+00 


.0 




.18 


.825378+00 


.0 




.20 


.809020+00 


.0 




.30 


.734559+00 


.0 




.40 


.670788+00 


.0 




.50 


.615690+00 


.0 




.60 


.567805+00 


.0 




.70 


.525930+00 


.0 




.80 


.489101+00 


.0 




.90 


.456532+00 


.0 




1 .00 


.427584+00 


.0 




1 .10 


.401730+00 


.0 




1 .20 


.378537+00 


.0 




1 .30 


.357643+00 


.0 


Q\ 


1 .40 


,33e744+00 


.0 


00 


1 .50 


.321585+00 


.0 


1 .60 


.305953+00 


.0 




1 .70 


.291663+00 


.0 




1 .80 


.278560+00 


.0 




1 .90 


.266509+00 


.0 




2 .00 


.255396+00 


.0 




2.10 


.245119+00 


.0 




2.20 


.235593+00 


.0 




2.30 


.226742+00 


.0 




2 .40 


.218499+00 


.0 




2.50 


.210806+00 


.0 




2.60 


.203613+00 


.0 




2 .70 


.196874+00 


.0 




2.80 


.190549+00 


.0 




2.90 


. 184602 + 00 


.0 




3.00 


.179001 +00 


.0 




3.50 


. 155294+00 


.0 




4 .00 


.136999+00 


.0 




4.50 


.122485+00 


.0 




5.00 


.1 10705+00 


.0 




5 .50 


.100962+00 


.0 




6.00 


.927766-01 


.0 




6.50 


.858057-01 


.0 




7 .00 


.798001-01 


.0 




7.50 


.745737-01 


.0 




8.00 


.655e52-01 


.0 




8.50 


.659251-01 


.0 




9.00 


.623077-01 


.0 



15° 


30° 




37.5 




45° 


100000+01 


.000000 


.100000+01 


000000 


.100000+01 


000000 


.100000+01 


.000000 


978544+00 - 


-.564511-02 


.980656+00 -. 


109433-01 


.982202+00 -. 


133574-01 


.984046+00 


-. 155619-01 


957755+00 - 


-. 109148-01 


.961711+00 -. 


212290-01 


.964623+00 -. 


259749-01 


.96811 7+00 


-. 303494-01 


937610+00 - 


-. 158323-01 


.943165+00 -. 


308907-01 


.947276+00 -. 


378843-01 


.952236+00 


-.443878-01 


918086+00 - 


-.204192-01 


.925014+00 -. 


399605-01 


.930171+00 -. 


491167-0 1 


.936422+00 


-.577025-01 


899159+00 ■ 


-.246960-01 


.907257+00 -. 


484692-01 


.913315+00 -. 


597018-0 1 


.920696+00 


-.703184-01 


880810+00 - 


-.286817-01 


.889890+00 -. 


564458-01 


.896717+00 -. 


696688-0 1 


.905075+00 


-.822605-01 


863016+00 - 


-.323942-01 


.872908+00 -. 


639183-01 


.880383+00 -. 


790455-01 


.889575+00 


-.935532-01 


845758+00 - 


-.358502-01 


.856308+00 -. 


709132-01 


.864318+00 -. 


878590-01 


.874212+00 


-. 104221+00 


.829017+00 


-.390656-01 


.840085+00 -. 


774557-01 


.848526+00 -. 


961353-01 


.858998+00 


-. 114286 + 00 


.812775+00 


-.420551-01 


.824233+00 -. 


835698-01 


.833010+00 -. 


103900+00 


.843946+00 


-. 123774+00 


.738441 +00 


-.540691-01 


.750340+00 -. 


108476+00 


.759640+00 -. 


135855+00 


.771453+00 


-. 163319+00 


.674177+00 


-.621855-01 


.684812+00 -. 


125660+00 


.693269+00 -. 


158259+00 


.704192+00 


-. 191604+00 


.618355+00 


-.674592-01 


.626815+00 -. 


137032+00 


.633647+00 -. 


173298+00 


.642609+00 


-. 210932+00 


.569635+00 


-.706620-01 


.575515+00 -. 


144060+00 


.580342+00 -. 


182718+00 


.586778+00 


-.223250+00 


.52691 1+00 


-.723596-01 


.530125+00 -. 


147863+00 


.532829+00 -. 


187900+00 


.536524+00 


-. 230166 + 00 


.489272+00 


-.729681-01 


.489917+00 -. 


149296+00 


.490547+00 -. 


189926+00 


.491521+00 


-.232989+00 


.455962+00 


-.727939-01 


.454239+00 -. 


149005+00 


.452940+00 -. 


189631+00 


.451358+00 


-.232766+00 


.426354+00 


-.720633-01 


.422513+00 -. 


147478+00 


.419480+00 -. 


187662+00 


.415588+00 


-.230320+00 


.399923+00 


-.709434-01 


.394231+00 -. 


145081+00 


.389679+00 -. 


184510+00 


.363760+00 


-.226295+00 


. 376233+00 


-.695576-01 


.368549+00 -« 


142090+00 


.363094+00 -. 


180548+00 


.355440+00 


-.221186+00 


.354916+00 


-.679971-01 


.346284+00 -, 


138710+00 


.339332+00 -. 


176056+00 


.320223+00 


-.21 5369+00 


.335662+00 


-.663291-01 


.325505+00 -« 


135093+00 


.318044+00 -• 


171 244+00 


.307738+00 


-.209128+00 


.318209+00 


-.646029-01 


.307525+00 -. 


131350+00 


.298923+00 -« 


166 267+00 


.287653+00 


-.202672+00 


.302334+00 


-.628545-01 


.290896+00 -. 


127564+00 


.281701+00 -. 


161 236+00 


.269672+00 


-. 196156+00 


.287848+00 


-.611100-01 


.275808+00 -« 


. 123794 + 00 


.266147+00 -. 


156235+00 


.253535+00 


-. 189690 + 00 


.274589+00 


-.593879-01 


.262075+00 - 


.120081 +00 


.252058+00 - 


151319+00 


.23901 3+00 


-. 183351+00 


.262415+00 


-.577015-01 


.249540+00 - 


. 116455+00 


.239260+00 - 


146528+00 


.225508+00 


-. 177191+00 


.251208+00 


-.560596-01 


.238065+00 - 


. 1 12935 + 00 


.227600+00 - 


.141 e89+00 


.214048+00 


-. 171246 + 00 


.240863+00 


-.544681-01 


.227532+00 - 


. 109534+00 


.216948+00 - 


.1 37417+00 


.203282+00 


-. 165535 + 00 


.231289+00 


-.529306-01 


.217837+00 - 


. 106258+00 


.207188+00 - 


133122+00 


. 193480 + 00 


-. 160068 + 00 


.222408+00 


-.514490-01 


.208891+00 - 


. 1031 12 + 00 


.198222+00 - 


.129008+00 


•ie4530+00 


-. 154849 + 00 


.214150+00 


-.500239-01 


.200615+00 - 


. 100095 + 00 


.189963+00 - 


.125073+00 


.176334+00 


-. 149876 + 00 


.206456+00 


-.486549-01 


.192541+00 - 


.972068-01 


.182336+00 - 


.121315+00 


.168808+00 


-. 145143 + 00 


.199271+00 


-.473413-01 


.185809+00 - 


.944437-01 


.175275+00 - 


•1 17730+00 


.161878+00 


-. 140642 + 00 


.152549+00 


-.460815-01 


.179166+00 - 


.918020-01 


.168723+00 - 


.1 14311+00 


.155479+00 


-. 136364 + 00 


. 186248+00 


-.448740-01 


.1 72967+00 - 


.892774-01 


.162630+00 - 


. 1 1 1051+00 


. 149556 + 00 


-. 132298 + 00 


.180331+00 


-.437168-01 


.167168+00 - 


.868650-01 


.156950+00 - 


. 107943 + 00 


.144061+00 


-. 128434 + 00 


.174766+00 


-.426079-01 


.161735+00 - 


.845599-01 


.151645+00 - 


. 1 04980 + 00 


. 138950+00 


-. 124761+00 


.151281 +00 


-.377161-01 


.139036+00 - 


.744669-01 


.129658+00 - 


.920841-01 


.1 17989 + 00 


-. 108897 + 00 


.133241+00 


-.337346-01 


.121846+00 - 


.663430-01 


.113196+00 - 


.817946-01 


.102526+00 


-.963790-01 


.1 18978+00 


-.304579-01 


.108404+00 - 


.597166-01 


.100432+00 - 


.734593-01 


.906646-01 


-.863253-01 


. 107435+00 


-.277278-01 


.976147-01 - 


.542350-01 


.902530-01 - 


.666012-01 


.812804-01 


-. 781076-01 


.979088-01 


-.254257-01 


.887699-01 - 


.496392-01 


.819492-01 - 


.608758-01 


.736702-01 


-. 712820-01 


.899189-01 


-.234628-01 


.813899-01 - 


.457386-01 


.750469-01 - 


.560327-01 


.673730-01 


-. 655313-01 


.831244-01 


-.217720-01 


.751399-01 - 


.42391 1-01 


.692192-01 - 


.518878-01 


.620749-01 


-.606252-01 


.772774-01 


-.203020-01 


.697797-01 - 


.394900-01 


.642331-01 - 


.483034-01 


.575547-01 


-.563934-01 


.721940-01 


-.190134-01 


.651322-01 - 


.369532-01 


.599186-01 - 


.451749-01 


.536520-01 


-.527075-01 


.677345-01 


-.178753-01 


.610646-01 - 


.347176-01 


.561483-01 - 


.424219-0 1 


.502479-01 


-.494694-01 


.637913-01 


-. 168633-01 


.574747-01 - 


. 327332-01 


.528254-01 - 


.399814-01 


.472522-01 


-.466031-01 


.602800-01 


-.159579-01 


.542832-01 - 


. 309606-01 


.498745-01 - 


.378036-01 


.445553-01 


-.440484-01 



EXP(Z**2J+ERFC(Z) 



ON 

CO 





RHOXTHETA 


50 




60 


70 


80 


90 






.00 


.1000004-01 


.000000 


.100000+01 


.000000 


. 100000+01 


.000000 


•100000+01 


.000000 


.100000+01 , 


000000 




.02 


.985430*00 


-.168969-01 


.988522+00 


-. 1 91978-01 


.991980+00 


-.209465-01 


.995708*00 


-.220828-01 


.999600+00 -« 


225616-01 




.04 


.970750+00 


-.330243-01 


.976680+00 


-.377036-01 


.983379+00 


-.413619-01 


.990684+00 


-.438614-01 


•998401+00 -< 


450871-01 




.06 


.955991+00 


-.484013-01 


.964508+00 


-.555200-01 


.974228+00 


-.612308-01 


.984947+00 


-.653065-01 


.996406+00 -. 


675405-0 1 




.00 


.941179+00 


-.630470-01 


.952039+00 


-.726507-01 


.964559+00 


-.805399-01 


.978518+00 


-.863902-01 


.993620+00 -. 


898862-01 




.10 


.926338+00 


-.76<S8l3-0l 


.939307+00 


-.891010-01 


.954404+00 


-.992775-01 


.971421+00 


-. 107C86+00 


.990050+00 -, 


112089+00 




.12 


.91 1494 + 00 


-.902242-01 


.926342+00 


-.104877+00 


.943794+00 


-.117434+00 


.963679+00 


-.127369+00 


.985703+00 -, 


134113+00 




. 14 


.896667+00 


-.102796+00 


.913175+00 


-.1 19987+00 


.932761+00 


-.135001+00 


.955317+00 


-. 147215+00 


.980591+00 -, 


155925+00 




.16 


.881678+00 


-.114717+00 


.899834+00 


-.134439+00 


.921336+00 


-. 151972+00 


.946360+00 


-.166602+00 


.974725+00 - 


.177491+00 




.18 


.867146+00 


-.126007+00 


.886348+00 


-. 148243+00 


.909551+00 


-• 168342 + 00 


.936834+00 


-.185511+00 


.968119+00 - 


.198777+00 




.20 


.852490+00 


-.136687+00 


.872742+00 


-.161410+00 


•897434+00 


-.184108+00 


.926768+00 


-.203921+00 


.960789+00 - 


.219753+00 




.30 


.760658+00 


-.181638+00 


.803774+00 


-.218107+00 


.632921+00 


-.253864+00 


.869300+00 


-.287964+00 


.913931+00 - 


.318916+00 




.40 


.713036+00 


-.214294+00 


.735131+00 


-.260834+00 


.764351+00 


-.308806+00 


.802581+00 


-.357672+00 


.852144+00 - 


406153+00 




.50 


.649979+00 


-.236926+00 


.668830+00 


-.291462+00 


.694712+00 


-.349876+00 


.730181 +00 


-.412479+00 


.778801+00 - 


478925+00 




• 60 


.592156+00 


-.251548+00 


.606260+00 


-.3 11920+00 


.626410+00 


-.378484+00 


.655426+00 


-.452678+00 


.697676+00 - 


•535713+00 




.70 


.539686+00 


-.259892+00 


.548281 +00 


-.324060+00 


.561255+00 


-.396308+00 


.581214+00 


-.479263+00 


.612626+00 - 


.576042+00 




.80 


.492447+00 


-.263412+00 


.495329+00 


-.329575+00 


.500496+00 


-.405125+00 


.509895+00 


-.493727+00 


.527292+00 - 


.600412+00 




.90 


.450161+00 


-.263301+00 


.447520+00 


-.329960+00 


.444890+00 


-.406682+00 


.443222+00 


-. 497872+00 


.444858+00 - 


.610142+00 


I 


.00 


.412460+00 


-.260519+00 


.404741+00 


-.326488+00 


.394786+00 


-.402602+00 


.382368+00 


-.493623+00 


.367879+00 - 


.607158+00 


I 


.10 


.378934+00 


-.255830+00 


.366725+00 


-.320215+00 


.350217+00 


-. 394324+00 


.327976+00 


-.482883+00 


.298197+00 - 


.593761+00 


1 


.20 


.349161+00 


-.249826+00 


.3331 1 1+00 


-.311994+00 


.310985+00 


-.383080+00 


.280243+00 


-.467415+00 


.236928+00 - 


.572397+00 


1 


.30 


.322730+00 


-.242965+00 


•303494+00 


-.302503+00 


.276740+00 


-.369883+00 


.239021+00 


-.448772+00 


.184520+00 - 


.545456+00 


1 


• 40 


.299256+00 


-.235593+00 


.277452+00 


-.292262+00 


.247038+00 


-. 355543+00 


.203909+00 


-.428255+00 


.140858+00 - 


.5151 13+00 


1 


.SO 


.278384+00 


-.227969+00 


.254575+00 


-.281666+00 


.221396+00 


-.340682+00 


.174350+00 


-.406909+00 


.105399+00 - 


.483227+00 


1 


.60 


.259793+00 


-.220280+00 


.234476+00 


-.271004+00 


.199324+00 


-.325767+00 


. 149701+00 


-.385528+00 


.773047-01 - 


.451284+00 


1 


.70 


.243198+00 


-.212662+00 


.216798+00 


-.260483+00 


. 180352+00 


-.311128+00 


.129300+00 


-.364688+00 


.555762-01 - 


.420388+00 


1 


.80 


.226349+00 


-.205210+00 


.201224+00 


-.250247+00 


.164045+00 


-.296992+00 


.1 12501+00 


-.344773+00 


.391639-01 - 


.391291+00 


1 


.90 


.215025+00 


-.197986+00 


.167470+00 


-.240389+00 


• 150014 + 00 


-.283502+00 


.987098-01 


-.326022+00 


.270518-01 - 


.364437+00 


2 


.00 


.203034+00 


-. 191033 + 00 


.175289+00 


-.230965+00 


.137912+00 


-.270739+00 


.873971-01 


-.308553+00 


.183156-01 - 


.340026+00 


2 


.10 


.192211+00 


-.184372+00 


.164467+00 


-.222006+00 


.127441+00 


-.258737+00 


.781048-01 


-.292400+00 


.121552-01 - 


■ 318073 + 00 


2 


.20 


. 18241 1+00 


-.178015+00 


.154819+00 


-.213521+00 


. 118345+00 


-.247499+00 


.704459-0 1 


-.277541+00 


.790705-02 - 


. 298468+00 


2 


.30 


.173510+00 


-.171964+00 


.146185+00 


-.205507+00 


.1 10409+00 


-.237007+00 


.641001-01 


-.263913+00 


.504176-02 - 


.281026+00 


2 


.40 


.165400+00 


-.166214+00 


.138430+00 


-.197951+00 


.1 03450+00 


-. 227227+00 


.588065-01 


-.251431+00 


.315111-02 - 


.265522+00 


2 


.50 


.157988+00 


-.160758+00 


.131439+00 


-.190836+00 


.973 162-0 I 


-.2181 18+00 


.543547-01 


-.240002+00 


.193045-02 - 


.251723+00 


2 


.60 


. 151 195+00 


-.155584+00 


. 125111+00 


-. 184138+00 


.918812-01 


-.209636+00 


.505775-01 


-.229526+00 


.115923-02 - 


.239403+00 


2 


.70 


. 144950 + 00 


-.150680+00 


•1 19364+00 


-.177833+00 


.870395-01 


-.201733+00 


.473424-01 


-.219915+00 


•682328-03 - 


.228355+00 


2 


.80 


.139193+00 


-.146031+00 


.114124+00 


-. 171899+00 


.827037-01 


-. 194364 + 00 


.445450-01 


-.211075+00 


.393669-03 - 


. 218399+00 


2 


.90 


.133871+00 


-.141624+00 


.109330+00 


-.166309+00 


.788013-01 


-. 187486 + 00 


.421033-01 


-.202926+00 


.222630-03 - 


.209377+00 


3 


.00 


.128940+00 


-.137445+00 


.104930+00 


-.161041+00 


.752717-01 


-. 181058+00 


.3995 30-01 


-.195396+00 


.123410-03 - 


.201 157+00 


3 


.50 


.108892+00 


-.1 19509+00 


.874505-01 


-.138780+00 


.617C8S-01 


-.154443+00 


.321214-01 


-. 165001 + 00 


.478512-05 - 


. 168830 + 00 


4 


.00 


.542797-01 


-.105481+00 


.750909-01 


-. 121741+00 


.525024-01 


-. 134607 + 00 


.270928-01 


-.143006+00 


.112535-06 - 


. 145954 + 00 


4. 


50 


.831655-01 


-.942896-01 


.658787-01 


-. 108356+00 


.45804 8-01 


-. 119295+00 


.235275-01 


-.126308+00 


.160523-08 - 


. 128735+00 


5, 


00 


.744259-01 


-.851875-01 


.587346-01 


-.975919-01 


.406890-01 


-. 107125+00 


.208412-01 


-.113170+00 


.138879-10 - 


. 115246 + 00 


5« 


50 


.673700-01 


-.776560-01 


.530232-01 


-.887581-01 


.3 66 39 9-01 


-.972205-01 


.187325-01 


-. 102546 + 00 


.728772-13 - 


. 104367+00 


6. 


00 


.615511-01 


-.713293-01 


.483467-01 


-.813834-01 


.333477-01 


-.890010-01 


.170271-01 


-.937703-01 


.231952-15 - 


.953962-01 


6. 


50 


.566679-01 


-.659441-01 


.444434-01 


-.751362-01 


.306137-01 


-.820692-01 


.156161-01 


-.863935-01 


.447773-18 - 


.878644-01 


7 . 


00 


.525101-01 


-.613075-01 


.41 1335-01 


-.697775-01 


.283042-01 


-.761436-01 


.144276-01 


-.801029-01 


.524289-21 - 


.814475-01 


7. 


50 


.489261-01 


-.572750-01 


.382898-01 


-.651309-01 


.263257-01 


-.710192-01 


.1341 15-01 


-.746733-01 


.372336-24 - 


.759126-01 


8. 


00 


.458040-01 


-.537368-01 


.358190-01 


-.61 0637-01 


.246107-01 


-. 665434-01 


.125322-01 


-.699379-01 


.160381-27 - 


•7108B1-01 


8 


50 


.430594-01 


-.506078-01 


.336516-01 


-.574741-01 


.2 3 1091-01 


-.626001-01 


.1 17632-01 


-.657709-01 


.419009-31 - 


.668445-01 


9 


,00 


.406273-01 


-.478214-01 


.317343-01 


-.542828-01 


.217828-01 


-.590993-01 


.110848-01 


-.620752-01 


.66396e-3£ - 


.630821-01 



TABLE 3 



EXP(Z**2»*ERFC(-Z) 



RHOMHETA 



On 
00 



.00 


.100000+01 


.0 


.02 


.102297+01 


.0 


.04 


.104678+01 


.0 


.06 


.107147+01 


.0 


.08 


.109708+01 


.0 


.10 


.1123C4+01 


.0 


.12 


.115122+01 


.0 


.14 


.117985+01 


.0 


.16 


.120958+01 


.0 


.18 


.124048+01 


.0 


.20 


.127260+01 


.0 


.30 


.145375+01 


.0 


.40 


. 167623+01 


.0 


.50 


.195236+01 


.0 


.60 


.229885+01 


.0 


.70 


.273870+01 


.0 


.80 


.330386+01 


.0 


.90 


.402928+01 


.0 


1 .00 


.500898+01 


.0 


1 .10 


.630524+01 


.0 


1 .20 


.806285+01 


.0 


1 .30 


.104813+02 


.0 


1 .40 


.138599+02 


.0 


1 .50 


. 186539 + 02 


.0 


1 .60 


.255657+02 


.0 


1 .70 


.356950+02 


.0 


1 .80 


.507889+02 


.0 


1 .90 


.736656+02 


.0 


2.00 


.108941+03 


.0 


2.10 


.164294+03 


.0 


2.20 


.252703*03 


.0 


2.30 


.396460+03 


.0 


2.40 


.634478+03 


.0 


2.50 


.103581+04 


.0 


2.60 


.172508+04 


.0 


2.70 


.293094+04 


.0 


2.80 


.508022*04 


.0 


2.90 


.898334+04 


.0 


3.00 


.162060+05 


.0 


3.50 


.417562+06 


.0 


4.00 


.177722*08 


.0 


4 .50 


.124593+10 


.0 


5.00 


.144010+12 


.0 


5.50 


.274434+14 


.0 


6.00 


.862246+16 


.0 


6.50 


.446655+19 


.0 


7.00 


.381469+22 


.0 


7.50 


.537149+25 


.0 


8 .00 


.124703+29 


.0 


8.50 


.477316+32 


.0 


9.00 


.301219+36 


.0 



15 


30 


37.5 


45 


.100000+01 


.000000 


.100000+01 


.000000 


.100000+01 


000000 


100000+01 


.000000 


.102215+01 


.604525-02 


.101974+0 1 


.116363-01 


.101801+01 


141302-01 


101595+01 


. 163619-01 


. 104502+01 


.125170-01 


.103989+01 


.240025-01 


.103620+01 . 


290671-01 • 


103188+01 


.335494-01 


.106863+01 


.194435-01 


.106043+01 


.371373-01 


.105458+01 . 


448455-01 . 


104775+01 


.515878-01 


.109302+01 


.268548-01 


.108137+01 


.510811-01 


.107311+01 


615C09-01 . 


106354+01 


.705024-01 


.111821+01 


.347829-01 


.1 10269+01 


.658763-01 


.109177+01 . 


790701-01 . 


107920+01 


.903181-01 


.114424+01 


.432622-01 


.1 12441 + 01 


.815669-01 


.111056+01 . 


975904-01 . 


109472+01 


. 111059 + 00 


.117112+01 


.523294-01 


.1 14650+01 


.981992-01 


.112943+01 • 


117100+00 , 


1 1 1004+01 


. 132751 + 00 


.1 19891+01 


.620234-01 


.1 16896+01 


. 115821+00 


.114836+01 . 


137638+00 • 


112513+01 


. 155415+00 


. 122763+01 


.723862-01 


.1 19178+01 


. 134483 + 00 


.116733+01 . 


159244+00 . 


113995+01 


. 179075+00 


.125731+01 


.834622-01 


.121495+01 


. 154237+00 


.118629+01 


181958+00 • 


115445+01 


.203753+00 


. 142149+01 


.151332*00 


.1 23536+01 


.271371*00 


.127976+01 


313595+00 . 


122045+01 


.343076+00 


.161572+01 


.245769+00 


.146100+01 


.424910+00 


.136645+01 . 


479143+00 • 


127026+01 


.510240+00 


.184573+01 


.377084+00 


.158657+01 


.623876+00 


.143813+01 • 


683551+00 . 


129522+01 


.705740+00 


.211791+01 


.559713+00 


. 170349+01 


.878536+00 


.148357+01 . 


930808+00 


128502+01 


.927798+00 


.243899+01 


.813903*00 


.179848+01 


.119994+01 


.148800+01 • 


122283+01 . 


122814+01 


. 117142 + 01 


.281529+01 


.1 16807+01 


.185201+01 


. 159889+01 


.143293+01 • 


155787+01 . 


111267+01 


. 142738+01 


.325121+01 


.166206+01 


.183636+0 1 


.208415+01 


.129634+01 . 


192846+01 


927639+00 


. 168134 + 01 


.374645+01 


.235168+01 


.171377+01 


.265934+01 


.105379+01 • 


231879+01 • 


665017+00 


. 191326+01 


.429102+01 


.331474*01 


.143482+01 


.331817+01 


.680682+00 . 


270192+01 


32227e+00 


. 209753+01 


.485657+01 


.465908+01 


.938062+00 


.403754+01 


.156290+00 • 


303701+01 


945930-01 


.220410+01 


.538160+01 


.653258+01 


.151541+00 


.476793+01 


-.530066+00 


326754+01 -i 


568066+00 


.2201 18+01 


.574679+01 


.913501+01 


-.998ei6+00 


.542135+01 


-.137052+01 


•332164+01 -. 


106664+01 


.205955+01 


.573426+01 


.127298*02 


-.257574+01 


.585743+01 


-.232809+01 


311621+01 -. 


154400+01 


.175882+01 


.496175+01 


. 176520+02 


-.462254+01 


.587039+01 


-.332543+01 . 


256681+01 -. 


194085+01 


. 129487+01 


.277771+01 


.243017+02 


-.708673+0 1 


.518191+01 


-.423536+01 


160541+01 -• 


219057+01 


.687583+00 


-.190182+01 


.331040+02 


-.98042 3+01 


.344908+01 


-.487786+01 i 


206798+00 -, 


222534+01 


-. 131466-01 


-.108404+02 


.443947+02 


-.1 24081+02 


.301777+00 


-.502963+01 - 


.157716+01 -, 


201049+01 


-.725740+00 


-.268411+02 


.581561*02 


-. 142543+02 


-.457055+01 


-.445375+01 - 


.358054+01 -i 


152134+01 


-. 134236 + 01 


-.542404+02 


.734660+02 


-.143607+02 


-. 112625 + 02 


-.295578+01 - 


.5494 15+01 -. 


798885+00 


-. 174372+01 


-.995202+02 


.874162+02 


-.1 14096+02 


-.194032*02 


-.468263+00 -• 


686146+01 i 


61 0500-01 


-. 182367+01 


-.171913+03 


.930880+02 


-.389141+01 


-.278220+02 


.284504+01 - 


.712231+01 


.907518*00 


-. 152069 + 01 


-.283609+03 


.759226+02 


.950566+01 


-.341808+02 


.649023+01 - 


.572745+01 


155613+01 


-. 849408+00 


-.448599+03 


.748935+01 


.291411+02 


-.347103+02 


.959591+01 - 


.233541+01 


.183009+01 


.787847-01 


-.677955+03 


-.164667+03 


.532363+02 


-.243320+02 


.109818+02 


.292473+01 < 


161505+01 


.105855+01 


-.967019+03 


-.532419+03 


.763516+02 


.239913+01 


.941085+01 


.919036+01 i 


913632+00 


. 182663+01 


-.126560+04 


-.124780+04 


.879744+02 


.489857+02 


.405857+01 


.147289+02 - 


.121594+00 


.213210+01 


-.141486+04 


-.254465+04 


.722456+02 


.112885+03 


-.4851 16+01 


.171054+02 - 


.119568+01 


. 182716+01 


-.102319+04 


-.474402+04 


.105895+02 


.179798+03 


-.154399+02 


.138266+02 -• 


196121+01 


.948998+00 


.759689+05 


-.1 27565+05 


-.345002+03 


-.846681+03 


.352504+02 - 


.3181 18+02 


.178275+01 


-.513342+00 


-.303148+06 


.206132*07 


.165184+04 


.572854+04 


-.121850+03 


.315751+02 - 


.201785+01 


-.479428+00 


-.632047+08 


-.532606+08 


.127488+05 


-.482630+05 


.286262+03 


•246430+C3 


.248394+00 


.205738+01 


.504457+10 


-.335304+09 


-.505861+06 


. 179234 + 06 


.714413+03 - 


.107588+04 


.190113+01 


-. 186596+00 


-.398070+12 


.262492+12 


.359254+07 


.647923+07 


-.294447+04 - 


,407345+04 


714128+00 


-. 176703 + 01 


.457872+14 


-.520744+14 


.127587+09 


-. 310891+08 


-.217456+05 - 


.476709+04 - 


.323300+00 


-.191803+01 


-.100716+17 


.1 18452 + 17 


.133025+10 


-.267648*10 


-.112172+06 


.340422+04 - 


.383674+00 


-. 191335 + 01 


.433489+19 


-.317887+19 


.207275+10 


-.873217+11 


-.630212+06 - 


. 132005+06 


.543630+00 


-. 185111+01 


-.283374+22 


.426347+21 


.630673+1 1 


-.327704+13 


-.252634+07 - 


.3361 14 + 07 


.185781+01 


-.535787+00 


.196497+25 


.129886+25 


.683728+14 


-. 142358+15 


.165524+08 - 


.265076+08 


.733467*00 


. 188952 + 01 


-.989764+25 


-.298524+28 


.943832+16 


-.252614+16 


.206701+09 


.164806+09 - 


.204721+01 


.598650-01 


-.549838+31 


.194934+31 


.397576+18 


.666613+18 


-.243169+10 


.751742+09 


.150878+01 


-. 121573 + 01 



EXPf 2**2l*ERFCf-Z» 





RMOVTMETA 


50 


60 


70° 


80° 


90° 




.00 


. 100000 fOl 


.000000 


.100000+01 


.000000 


.100000+01 


.000000 


. 100000+01 


.000000 


.100000+01 


.000000 




.02 


.101443*0 1 


.176847-01 


. 101 108+01 


. 19B904-01 


.100741+01 


.214606-01 


.100354+01 


.223563-01 


.999600+00 


.225616-01 




.04 


.102869*01 


.36 1748-01 


.102172+01 


.404727-01 


.101417+01 


.434163-01 


.100631+01 


.449542-0 1 


.598401+00 


.450871-01 




.06 


.104275*01 


.554874-01 


.103189+01 


.617442-01 


. 102026 + 01 


.658462-01 


. 100830 + 01 


.677607-01 


.996406+00 


.675405-01 




.08 


. 105656*0 1 


.736385-01 


.104154+01 


.837004-01 


.1 02564+01 


.887273-01 


.100949+01 


.907418-01 


.993620+00 


.898862-01 




.10 


. 107009*01 


.966430-01 


.105064+01 


.106335+00 


.103029+01 


.112035+00 


.100986+01 


. 1 13862 + 00 


.990050+00 


. 112089 + 00 




.12 


.108331*01 


. 118515*00 


.105916+01 


.129639+00 


. 103418 + 1 


. 135743+00 


. 100942+01 


. 137086 + 00 
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